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We study correlated two-level quantum dots, coupled in effective 1-channel fashion to metallic leads; with 
electron interactions including on-level and inter-level Coulomb repulsions, as well as the inter-orbital Hund's 
rule exchange favoring the spin-1 state in the relevant sector of the free dot. For arbitrary dot occupancy, the 
underlying phases, quantum phase transitions (QPTs), thermodynamics, single-particle dynamics and electronic 
transport properties are considered; and direct comparison is made to conductance experiments on lateral quan- 
tum dots. Two distinct phases arise generically, one characterised by a normal Fermi liquid fixed point (FP), the 
other by an underscreened (USC) spin-1 FP. Associated QPTs, which occur in general in a mixed valent regime 
of non-integral dot charge, are found to consist of continuous lines of Kosterlitz-Thouless transitions, separated 
by first order level-crossing transitions at high symmetry points. A 'Friedel-Luttinger sum rule' is derived and, 
together with a deduced generalization of Luttinger's theorem to the USC phase (a singular Fermi liquid), is used 
to obtain a general result for the T = zero-bias conductance, expressed solely in terms of the dot occupancy 
and applicable to both phases. Relatedly, dynamical signatures of the QPT show two broad classes of behavior, 
corresponding to the collapse of either a Kondo resonance, or antiresonance, as the transition is approached from 
the Fermi liquid phase; the latter behavior being apparent in experimental differential conductance maps. The 
problem is studied using the numerical renormalization group method, combined with analytical arguments. 



I. INTRODUCTION 

The Kondo effect is one of the enduring paradigms of quan- 
tum many-body theory^. For most of its history it has been as- 
sociated with bulk condensed matter, notably transition metal 
impurites dissolved in clean metals, and certain heavy fermion 
rare earth compounds^. In recent years, however, the ad- 
vent of quantum dot systems - with the impressive control 
and tunability possible for 'artificial atoms' - has generated a 
strong resurgence of interest in Kondo and related physics in 
nanoscale devices (for reviews see e.g. JUS])- 

In odd-electron quantum dots the spin- 1/2 Kondo ef- 
fect arises. Manifest experimentally^ as a strong low- 
temperature enhancement of the zero bias conductance, indi- 
cating the formation of the local Kondo singlet below a char- 
acteristic Kondo temperature, the basic theoretical model here 
is of course the Anderson impurity model&: a single dot level, 
with a single on-level Coulomb interaction, tunnel coupled to 
non-interacting metallic leads. Moreover the Anderson model 
captures not only the Kondo regime - arising towards the cen- 
ter of the associated Coulomb blockade valley where the dot 
level is singly occupied - but also the mixed valent regimes 
of non-integral occupancy occurring towards the edges of the 
valley. As such, it encompasses essentially all the physics as- 
sociated with a single 'active' dot level. 

The situation is naturally more complex, and richer, if two 
active dot levels are integral to electronic transport. For exam- 
ple, higher dot spin states now become possible, in this case 
a 2-electron triplet stabilised by the inter-orbital Hund's rule 
exchang This state has been observed experimentally in 
even-electron dots, for both latera l 9 ! 10 ! 11 ' 12 and vertical^ de- 
vices (as well as in a single-molecule dot-). It too is manifest 
in a strong enhancement of the zero-bias conductance, indica- 
tiv e) 15 ! 16 of proximity to an underscreened spin-1 fixed point 17 
in which the spin-1 is quenched to an effective spin-^ on cou- 
pling to the leads. 



Much important theoretical work on the problem has en- 
sued; inclu ding both the 1-channel case (see e.g. Il6lll8lll9l 
l20H2lE l23l24ll ) where the single screening channel yields 
underscreened (USC) spin-1 as the stable low-temperature 
fixed point, and the 2-channel cas o 15 ' 16,23 - 25 ! 26 ' 27 where the 
spin-1 local moment is fully screened at the lowest temper- 
atures. Further, since the USC spin-1 fixed point is clearly 
distinct from that characteristic of a normal Fermi liquid - 
the USC phase being a 'singular Fermi liquid' 28 - quantum 
phase transitions from a normal Fermi liquid to the USC 
phase are expected, and found, to arise in the 1-channel case 
(with pristine transitions broadened into crossovers for 2- 
channel screening). This too has been studied quite exten- 
sivel y 1 8 i 1 V Q i 22 - 24 - 25 i 26 i 27 . However the large majority of pre- 
vious work on these 'singlet-triplet' transitions has focussed 
on a somewhat particular case - the middle of the 2-electron 
Coulomb blockade valley where, throughout both phases, the 
dot occupancy /charge remains close to 2; a situation we regard 
as unlikely to be applicable to a transition driven by tuning a 
gate voltage in the absence of a magnetic field (as in the exper- 
iments e.g. of ifllll ). where one instead expects the dot charge 
to vary continuously with gate voltage. A notable exception is 
the work of l22ll . in which low-temperature transport is con- 
sidered in a region separating two adjacent Coulomb blockade 
valleys with spins S = \ and S = 1 on the dot, and where the 
resultant quantum phase transition, driven by gate voltage and 
arising in the limit B — > 0+ of vanishing magnetic field, oc- 
curs in a mixed valent regime of non-integral dot charge. 

In view of the above our aim here is to consider a rather 
general model of a two-level quantum dot, coupled in a 1- 
channel fashion to metallic leads; to consider its underlying 
phases, thermodynamics, single-particle dynamics and asso- 
ciated low-temperature (T) electronic transport, for arbitrary 
dot charge - spanning as such the full range of possible be- 
havior; and ultimately to make tangible comparison to experi- 
ment^. The model itself is specified in sec.|II]and reflects the 
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natural complexity of a two-level dot, where in addition to the 
one-electron dot levels electron interactions include both on- 
level and inter-level Coulomb repulsions, together with inter- 
orbital spin-exchange. We study it using Wilson's numerical 
renormalization group (NRG) techniqu e 29 ' 30 ' 31 as the method 
of choice, employing the full density matrix formulation of 
the metho d 32 ' 33 (for a recent review see |34]); together where 
possible with analytical arguments. 

The intrinsic phases and associated thermodynamics are 
considered in sec. [HI] With Ei, £2 denoting the one-electron 
level energies, the general structure of the phase diagrams in 
the (£1 , £2)-plane is found to consist of a closed, continuous 
line of quantum phase transitions (QPTs) separating an USC 
spin-1 phase from a continuously connected normal Fermi 
liquid phase; although more complex topologies arise as the 
exchange coupling is driven weakly antiferromagnetic (sec. 
IIII B\ . leading ultimately to destruction of the USC phase. The 
transitions are found in general to be of Kosterlitz-Thouless 
type, except for particular lines of symmetry where first order, 
level-crossing transitions arise (sec. IIII At . 

Sec. [IV] focusses on the T = zero-bias conductance G c , 
and associated static phase shift 8. A 'Friedel-Luttinger sum 
rule' for 8 is derived, applicable to both the normal Fermi 
liquid and the USC spin-1 phases, and reducing to the usual 
Friedel sum rule-l^ in the Fermi liquid phase. Since the USC 
phase is a singular Fermi liquid 28 , and as such not perturba- 
tively connected to the non-interacting limit of the model, one 
does not expect Luttinger's (integral) theorem^ to apply. A 
generalization of it for the USC phase is however deduced, 
and its important consequences for the zero-bias conductance 
considered; leading to a simple result which, for both the nor- 
mal Fermi liquid and USC phases, gives G c in terms of the dot 
occupancy/charge (or, strictly, the 'excess impurity charge'-). 

Single-particle dynamics for both phases are detailed in sec. 
fVl In particular, dynamical signatures of the QPT on ap- 
proaching it from the normal Fermi liquid are found to fall into 
two broad classes, corresponding respectively to an 'on the 
spot' vanishing of either a Kondo resonance, or a Kondo an- 
tiresonance, in the single-particle spectrum; the spectral col- 
lapse in either case being associated with a vanishing Kondo 
scale Tg as the transition is approached, and in terms of which 
universal scaling of dynamics is found to occur. 

Finally, in sec.[VT]we make explicit comparison to the ex- 
periments of iflTll on a lateral dot; in which, on continuous 
tuning of a gate voltage at zero magnetic field, both the nor- 
mal spin-1 /2 Fermi liquid and the USC spin-1 phase are ob- 
served in adjacent Coulomb blockade valleys. Both the zero- 
bias conductance as a function of gate voltage, and (in this 
case inevitably approximate) differential conductance maps as 
a function of both gate and bias voltages, are compared to ex- 
periment; and the features observed related to the dynamics 
considered in sees. dIVIVl . We believe it fair to say that the 
underlying theory accounts rather well for experiment. 



II. MODEL 

Interacting quantum dots and other nanodevices are de- 
scribed generally by the dot H^, a pair of non-interacting 
leads Hl, and a tunnel coupling between the subsystems: 
H = Hp +Hl + Hj. We consider in this work a two-level 
interacting quantum dot of form: 

h d = Jl(e i + ^Un ( _ a )n ia + U'n l n 2 -J I1 n 1 -%2 (1) 

i.a 

Here h ia = d] a d ia where d] a creates a a (=t,l) spin elec- 
tron in level ; (= 1,2), h i = L CT ra !Cr is the total number oper- 
ator for level i, and s, is the local spin-operator with compo- 
nents sf = Y.a,a' dfaO^idj 1 and the Pauli matrices. The 
single-particle levels have energies £,-, the on-level Coulomb 
interaction (taken to be the same for both levels) is denoted 
by U, and the inter-level interaction by U'. Finally /h is the 
exchange coupling, taken in accordance with Hund's rule to 
be ferromagnetic (/h > 0, although we also comment in sec. 
IIIIBI on the weakly antiferromagnetic case). The states arising 
from He, itself will be discussed in sec. IIIBI below. 

The Hamiltonian for the two equivalent non-interacting 
leads (V = L,R) is given by H L = £ v L k , a £ k fl^ V(J a kV(J . Tun- 
nel coupling to the leads is described generally by Hf = 

Lv L,k,a V /v(4r a kv(T + h - c -) where v iv is the ^nel coupling 
matrix element between dot level i and lead v. We consider 
explicitly in this paper the case of an effective 1 -channel setup, 
in which the ratio V 2v /V lv = V 2 /Vj is independent of the lead 
index v; i.e. the tunnel couplings are of form \ iL = «V ( -, 
W jR = J3V ; (with a 2 + fi 2 = 1), as illustrated schematically 
in fig. Q] A simple canonical transformation to new lead 
orbitals may then be performed, c ka = tta^ + fia^^ and 
c ka = —j3a kL(J + aa kRa , such that solely the bonding combi- 
nation of lead states (c k(J ) couples to the dot: 

Ht = E V ; .(4c ka + h.c.) (2) 

(,k,(T 

We can thus drop the lead index v and consider one effective 
lead 

h l = EvU,. ( 3 ) 

hence the effective 1 -channel description illustrated in fig. [T] 
In practice we consider the standard casei of a symmetric, 
flat-band conduction band with half bandwidth D, i.e. the lead 
density of states (per conduction orbital) is 1 /(2D). 

It need hardly be added that the tunnel coupling pattern 
considered (fig.[U) is not the most general case, which would 
by contrast involve an irreducibly 2-channel description^, in 
general with strong channel anisotrop y 15 ' 16 . The richness of 
the physics arising in the case considered, with its associ- 
ated pristine QPTs, is nonetheless more than ample to justify 
its study (indeed for NRG calculations in practice, we focus 
largely on the case V 2 = Vj). It has moreover been argued 
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FIG. 1: Schematic of the two-level dot (levels 1,2) coupled to the V = 
L,R leads with the tunnel couplings V^v/^lv — V^/Vi independent of 
v; leading to the equivalent one-lead description indicated. 



A. Symmetries 

We will subsequently consider different phases of the dot- 
lead coupled system in the (£i,£2)-plane, for given values of 
the interaction parameters U, U' and 7h entering (eqn. (Q])). 
To this end it is economical to exploit symmetry. Rather than 
the bare levels E\ , E 2 it is often helpful to employ 



u' 



(11a) 



(see e.g. 1 1 9||2(Jtl ) that the 1 -channel case is generally appro- 
priate to lateral quantum dots, while a 2-channel model is ap- 
propriate for vertical dots. 

In considering equilibrium electronic transport per se, the 
central quantity is of course the zero-bias conductance, G C (T), 
across the L/R leads (fig. [TJ. An expression for it is readily 
obtained following Meir and Wingreen2I, and is: 

-dm 



Gc(T) = ^-G 
h 



dco 



da 



7i(Tn+r 2 2)D ee ((o) (4) 



Here /(ft)) = [e m l T + 1] 1 (kg = 1) is the Fermi function and 



(5) 



is the hybridization strength for level i (p is the lead den- 
sity of states). The dimensionless conductance prefactor 
Go = (2ap) 2 - or equivalent^ G = 4r L r fi /(r L + r R ) 2 with 
r v = Kp Y,i V iv - reflects the relative asymmetry in tunnel 
coupling to the L/R leads. It is naturally maximal, Go = 1, 
for symmetric coupling where (fig. [TJ a = 1/V2 = j5 (i.e. 
Fr = Fl). The key quantity determining the conductance eqn. 
©, which we analyse in detail in later sections, is the 'even- 
even' single-particle spectrum: D ee (a>) = — ^lmG ee ((o) in 
terms of the (retarded) Green function G ee ((o) («-> G ee (t) = 
— i6(t)({d ea (t),dl a })). The e-orbital creation operator is 
given generally by 

1 



dL 



F u d 



^22 dl 



VMT+T22 

in terms of the level creation operators, such that 

1 



2a 



G e e((0) = 



^22 f 



(6) 



(7) 



in terms of the corresponding propagators for the dot levels, 
Gjj((o) (ij £ {1,2}); and where 
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npV l W 2 = ^F n F 2 2 (8) 

is the inter-level hybridization strength. 

For the case V 2 = Vj all hybridization strengths coincide, 



•j 



V.=V 



(9) 



(and the e-e propagator then reduces simply to G ee ((o) = 
5 [Gn (ft)) + G22(ft>) +2Gi2(ft>)])- It is convenient in this case 
to specify the 'bare' parameters of Ho in terms of F, defining 



r' r' 
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.v = £2 



\u + u' 



(lib) 



Their significance arises from a particle-hole transformation 
(p-ht) of H = H D +H L + H T (eqns. ([TJ-O), namely^ 



die G 



ka 



k(j 



(12) 



H = H(x,y) transforms under the p-ht as H(x,y) — * 2(x+y) + 
n(—x, —y), and is hence invariant at the p-h symmetric point 
x = = y. Use of x,y thus specifies the level energies rel- 
ative to this point. All physical properties, thermodynamic 
and dynamic, have characteristic symmetries under the p-ht , 
which we exploit many times in the paper. For example the 
free energy F(x,y) = — TlnjTr e~£"w)} is, modulo an ir- 
relevant constant, equivalent to its p-ht counterpart (F(x,y) = 
2(x+y) + F(— x, — y)), whence e.g. phase boundaries (sees. 
(Ill BIIIH ) are invariant under inversion (x,y) — > (—x, —y)', and 
thus only y > x need in practice be considered. 

The second symmetry exploited is a '1-2' transformation, 
viz. the trivial canonical transformation 



( d 2cAc 



(13) 



under which the dot Hamiltonian Hjy(x,y) — > Ho(y,x). The 
same symmetry applies to the full H for V 2 = Vp whence 
e.g. F(x,y) = F(y,x) is invariant to reflection about the line 
y = x, and in consequence phase boundaries need overall be 
considered only fory >\x\. 



B. Ground state phases: overview 

It is first instructive to consider briefly the states of the iso- 
lated dot in the (x,y)-plane, as determined by the ground states 
of Ho (eqn. (UJ). We label the dot states as («i,«2) (with 
Hi = (hj ) the ground state charge for level /), with energies 
Er>(n\,n-i). For all /h > 0, the 2-electron dot ground state is 
the (1,1) spin triplet with energy Ej (= Eo{ 1 , 1 ) = £1 + £2 + 
U' — jJh), centred on the p-h symmetric point (x,y) = (0,0); 
as illustrated in fig. [2 a) for the representative case U = 20, 
U' = 7.5 and /h = 5. All other ground states, indicated in the 
figure, are either spin singlets or doublets. 

Considering in particular y > \x\ - phase boundaries be- 
ing invariant to inversion and reflection as above - the (1,1) 
triplet is bordered both by another 2-electron state, viz. the 
spin singlet (2,0), and by a 1 -electron spin doublet (1,0) (the 
bounding lines for which are given byy=x+U-U' + y H 
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FIG. 2: (a)-(c) show ground states of the isolated dot in the (x, v)-plane, for U = 20 and Jh =5 (in units of T = 1); with (a) U' = 7.5, (b) U' = 
and (c) U' = U. States are labelled as («i,«2), with n; = (hj) the charge on level ;'. Phase boundaries surrounding the 2-electron triplet state 
(«i ,«2) = (1, 1)) are indicated by solid lines, and all others by dotted lines. The dashed arrowed line in (a) shows the form of an experimental 
'trajectory' on application of a gate voltage, see text. For the same bare parameters as (a), Fig.[2{d) shows the phase diagram obtained via 
NRG for the lead-coupled dot system (with Vo = V\), as detailed in sec. Hill ff. It consists of a line of continuous quantum phase transitions 
(thick solid line) and two first order level-crossing transitions on the line y = x (shown as dots), separating a singular Fermi liquid phased 
characterised by an underscreened spin-1 fixed pointH (interior, 'USC') from a normal Fermi liquid ('FL') characterised in general by a frozen 
impurity FP— The hexagonal boundary of the (1, 1) triplet state for the isolated dot (fig.|2ja)) is also shown for comparison. 



and y = + jJh respectively). The dashed line in fig.|2a) 
shows a typical 'trajectory', y = x + Ae {i.e. £2 = £1 + Ae), 
expected from experiment on application of a gate voltage V g 
to the dot, with £\ y g and fixed level spacing Ae = £2 — £\ . 
Note that the 2-electron (1,1) triplet is thereby accessed from 
the 1 -electron state (1,0)^ (as relevant to comparison with 
experiment, sec. IV1I >. 

Fig. |2jb,c) show the isolated dot ground states arising for 
inter-level Coulomb repulsion U' = and U' = U, respec- 
tively. For y > \x\ in the former case, the (1,1) triplet state 
is bordered almost exclusively by the 1 -electron state (1,0), 
while in the latter case it is bordered almost exclusively by 
the 2-electron singlet (2,0). These two cases are of course 
extremes; and although aspects of the model have been con- 
sidered previously for the case U' = t / 18 i 19 ' 22 , we know of no 
compelling reason why the intra- and inter-level Coulomb re- 
pulsions should in general be near coincident for reasonably 
small dots (indeed we argue in sec.[VT]that comparison to the 
experiment of Hill is consistent with the contrary). 

The states of the dot per se are of course quite trivial. We 
now consider the full lead-coupled system, our aim here being 
to give simple qualitative arguments for the general form of 
the phase diagram in the (x,y)-plane. 

On coupling to the leads, the effective low-energy model 
deep in the spin-1 regime centred on (x,y) = (0,0), is nat- 
urally a 1 -channel spin-1 Kondo mode l 15 ! 16 (obtained for- 
mally by a Schrieffer- Wolff transformationi^ retaining only 
the triplet (1, l)-state of the dot itself, see also Appendix A). 
Its low-energy physics is well knownii half the spin-1 is 
screened by the conduction electrons, leading to a free spin- 
j with weak residual ferromagnetic coupling to the metallic 
lead, which results in in non-analytic (logarithmic) corrections 
to Fermi liquid behavior; the resultant state being classified as 
a singular Fermi liquid 2 ^. The associated low-energy fixed 
point (FP) is of course the underscreened spin-1 (USC) FP of 
Nozieres and Blandinii. 



By contrast, deep in the 1 -electron (1 , 0)-regime (fig. [2]), the 
effective low-energy model is obviously spin-^ Kondo, a nor- 
mal Fermi liquid with a fully quenched spin and a strong cou- 
pling (SC) low-energy F P 1 ! 29 ' 30 . Since the underlying stable 
FPs (USC and SC) associated with these two regimes are fun- 
damentally distinct, a quantum phase transition (QPT) some- 
where between the two must therefore occur—. 

But what of the other isolated dot states, encircling the spin- 
1 state as illustrated in fig. 0a-c) (and all of which as noted 
above are either spin singlets or doublets)? The salient point 
here is that, on coupling to the leads, all such give rise to 
Fermi liquid states: their stable low-energy FPs form a contin- 
uous line connecting the SC FP arising for the spin-| Kondo 
model to the generic case of the frozen impurity FP 3 ^ (as fol- 
lows from the original work of Krishnamurthy, Wilkins and 
Wilson 31 on the asymmetric single-level Anderson model). 
No phase transitions between these states can therefore oc- 
cur, the 'transitions' arising in the isolated dot limit (dotted 
lines in fig.|2a-c)) being replaced by continuous crossovers. 

In consequence, one expects the general structure of the 
phase diagram in the (x,y)-plane to consist of a closed, con- 
tinuous line of QPTs separating an USC spin-1 phase from a 
continuously connected normal Fermi liquid phase. This is 
indeed as found from detailed NRG analysis, as will be seen 
in the following sections. A typical resultant phase diagram is 
shown in fig. |2jd) (for the same bare parameters as fig. a), 
the phase boundary occurring close to the border of the (1,1) 
state of the isolated dot as one might expect). It consists of a 
line of continuous QPTs; together with two first order level- 
crossing QPTs on the line y = x (indicated by dots in fig.|2jd)), 
which are equivalent to each other under the p-h transforma- 
tion x — ► — x. 

The transitions will be discussed in detail below, but we 
add here that the occurrence of first order transitions along 
the y = x line (£2 = £1) is a general consequence of symme- 
try. As noted in sec. Ill Al for V2 = V\ the full H transforms 
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under the '1-2' transformation as H(x,y) — > H(y,x), and is 
hence invariant on the line y = x. Along that line all states 
of the entire system thus have definite parity under the '1-2' 
transformation, with the Hilbert space of H strictly separable 
into disjoint parity sectors. A level-crossing transition must 
thus occur when the global many-body ground state changes 
parity (further discussion of it will be given below). 



m. PHASES AND THERMODYNAMICS 

Dynamics and transport properties will be discussed in sec. 
|IV| ff, but we begin with thermodynamics; in particular the 
temperature (T) dependence of two standard quantities!^ 
which provide clear signatures of the various FPs reached 
under renormalization on decreasing the temperature/energy 
scale, namely the entropy Sj mp (r) and the uniform spin sus- 
ceptibility Ximp ( T ) = ( (S z ) 2 ) i mp / T (where S z refers to the spin 
of the entire system, and (£2)i mp = (£2) — (£2)o with (fl)o de- 
noting a thermal average in the absence of the dot). 

We also consider briefly the usual T = 'excess impurity 
charge' «j mp , viz. the difference in charge of the entire sys- 
tem with and without the dot present (£2 = N = Y.k a c k<r c k<T + 
J^n. in the above); and which in practice corresponds closely 
to the net dot charge, nj mp ~ (h\ see also sec. IIV A| 

Prosaic though n lmp is, we show later that it plays a key role 
in understanding the zero-bias conductance in both the USC 
and FL phases, and relatedly the 'Friedel-Luttinger sum rule' 
of sec. IIV Bl Under the p-h and 1-2 transformations of sec. 
ha] «i m p = rii m p(x,y) transforms respectively as: 
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Results shown are obtained using the full density ma- 
trix formulation 32 ^ of Wilson's non-perturbative NRG tech- 
niqu e) 29 ' 30,31 , employing a complete basis set of the Wilson 
chain; for a recent review see J3. Calculations are typically 
performed for an NRG discretization parameter A = 3, retain- 
ing the lowest 2000 states per iteration. We here consider ex- 
plicitly the case V% = V\ (sec. [TTJ, with the hybridization Y 
(eqn. (O) as the basic energy unit, choosing the lead band- 
width D/Y = 100 (3> 1, such that results are independent of 
D for all practical purposes). 

Fig. |3] shows the T/r-dependence of 5i mp (r) (top) and 
TXim P (T) (bottom), for fixed = 20, U' = 7 and 7 H = 2, tak- 
ing a vertical cut through the (jc,y)-phase diagram: the energy 
of level- 1 is fixed at l\ = -\tJ - U' = -17 (i.e. x = 0), and 
£2 (or equivalently y) is progressively decreased from deep in 
the FL phase, towards and through the transition, and down 
to the p-h symmetric point at the center of the USC phase; 
the transition occurring at 6^ = —6.536... (close to the value 
—6.5 expected from the isolated dot limit). 

In all cases the highest T behavior is naturally governed 
by the free orbital F P 3Q i 31 , with all 4 2 states of the two-level 
dot thermally accessible, hence S\ mp — In 16 (and T%\mp = 
2 x A). For case (a), £2 = +5 is sufficiently large that 
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FIG. 3: T-dependence of the entropy Sj mp (r) (top panel) and spin 
susceptibility r^j mp (r) (bottom); for fixed level energy Si =—jU — 
U' (i.e. x = 0) on progressively decreasing h from deep in the FL 
phase, across the QPT (at £2 C = —6.536..), and through to the p- 
h symmetric point atthe center of the USC phase. Shown for U = 
20, 0' = 7 (i.e. E\ = -17) and / H = 2, with: e 2 = +5 (a), -5 (b), 
-6 (c), -6.3 (d), -6.43 (e), -7 (f), -10 (g) and -17 (h). [Labels 
f) and g) are omitted for clarity from the top panel, but are easily 
identified from the bottom.] Inset, top panel: the corresponding 
T = impurity charge «; mp vs £ 2 . It changes continuously as the QPT 



is crossed at £2 C (dashed vertical line), at which point /; 



imp 



1.4. 



level-2 is in essence irrelevant (provided T/Y <C £2 — SiX 
the model thus reducing in effect to a single-level Anderson 
mode l 30 ! 3 1 . Hence, on decreasing T, 5i mp (r) first flows to- 
wards the spin- 1/2 local moment (LM) FP corresponding to 
Simp = In 2 (evident in this case as a relatively weak plateau 
at T/Y < 1, reflecting the modest minimum thermal excita- 
tion of - £d(2,0) -E d (1,0) = £i+f/ = 3Y). On further 
decreasing T, the system then flows to the stable frozen im- 
purity (FI) FP symptomatic of the Fermi liquid ground state, 
with vanishing entropy 5; mp (likewise T% lmv = 0). A Kondo 
scale Tk may be identified from the crossover between the 
marginally unstable LM FP and the stable FI FP (we define 
it in practice via Sj mp (7jr) =0.1). On further decreasing e%, 
cases (b-e) in fig. [3] the same essential behavior is found, the 
FIFP remaining the stable low-7^ FP. But the 5im P — ln2 (and 
TXimp — 4) LM plateau progressively lengthens and the asso- 
ciated Tk correspondingly diminishes, vanishing as the transi- 
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FIG. 4: Evolution of the Kondo scale in the FL phase £2 > &i c for the 
same parameters as fig. [5] log(T^/r) vs Eo — he- Inset: Tk vanishes 
exponentially as the QPT is approached, Tk °c exp(— 0-/\^2~ ^2cl) 
(with a ~ ^(1)), characteristic of a Kosterlitz-Thouless transition. 



tion is approached from the Fermi liquid side (fig. HJ. 

The behavior on the other side of the transition 62 < &i c 
(cases (f-h) in fig. |2]i is qualitatively distinct. Here the T = 
entropy is in all cases ln2 (with TXimp = \), characteristic 
of an unquenched doublet ground state. The stable FP is the 
spin- 1/2 LM FP - or equivalently the USC spin-1 FP— , there 
being no distinction between them as FPs per se. 

The QPT itself is of Kosterlitz-Thouless (KT) type. This is 
evident for example from NRG flows, which indicate no sep- 
arate, unstable critical FP, distinct from one of the stable FPs 
mentioned above. It is also evident in the behavior of the scale 
Tk, which as shown in fig. [4] (inset) vanishes exponentially in 
\&i — B2c\ (rather than as a power-law) as the QPT is ap- 
proached from the Femi liquid side; and by the absence of a 
low-energy scale in the USC phase which vanishes as the tran- 
sition is approached from that side. We add that the latter does 
not of course imply the inherent absence of a low-energy scale 
in the USC phase. For deep inside this phase (where «; mp ~ 2) 
the effective low-energy model is spin-1 Kondo, as evident 
e.g. in case (h) of fig.[3]from both the emergence of a near free 
spin-1 susceptibility with decreasing T (Ximp ^ S(S + 1)/3T 
with 5=1), and from the intermediate S lmp = ln3 plateau 
indicative of a spin-1 local moment FP; reached before the 
crossover to the stable USC FP with Sjmp = ln2, and from 
which a characteristic spin-1 Kondo scale 7J? = 1 may be identi- 
fied (in parallel to that above for the Fermi liquid Kondo scale 
Tk). But this scale plays no role in the QPT per se, and in con- 
trast to the approach from the Fermi liquid phase, there is no 
vanishing scale on approaching the QPT from the USC side. 

The behavior outlined above is not confined to the example 
illustrated: all continuous transitions are found to be of KT 
form. This is in fact to be expected. Hofstetter and Schoeller— 
have consider the model (with U' = U) in the regime where 
the dot is doubly occupied by electrons, i.e. «j mp ~ 2 - where 
from eqn. ( fT4l , n lmp = 2 arises by symmetry along the line 
y = —x in the phase-plane (or close enough to it, in practice). 
A KT transition is likewise found in this case-^, and by conti- 
nuity one would thus expect the same behavior to arise gener- 
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FIG. 5: Phase boundaries in the (x,y)-plane, separating the Fermi 
liquid phase (outer) from USC spin-1 (inner). Top: For U = 20 and 
In = 5, varying U' = (a), 7.5 (b) and 20 (c). First-order transitions 
on the line y = x are indicated by dots. Bottom: For U = 1 and U' = 
5.25, varying the exchange coupling Ju = 5, 1,0 and —0.5 (outside 
to inside respectively). Note the continued persistence of the USC 
spin-1 phase, even for weak antiferromagnetic exchange (see text). 



ally in the (x,y)-p\ane. 

We also note that the transition itself occurs generically in 
a mixed-valent regime of non-integral n lmp ; see e.g. fig. [3] 
(top, inset) where «; mp varies continuously as the transition 
is crossed, with «; mp ~ 1 .4 at the transition itself. This in turn 
means that even in a strongly correlated regime it is not in 
general possible to construct, via a Schrieffer-Wolff^ (SW) 
transformation from the original Anderson-like model, an ef- 
fective low-energy spin model applicable in the vicinity of the 
QPT. An exception to this is the vicinity of the line y = —x 
along which, as above, «ipm = 2 is guaranteed by symmetry. 
In this case, as shown in Ill9ll . a SW transformation retaining 
solely the 2-electron (1,1) triplet and (2,0) singlet states of 
the isolated dot yields an effective two-spin, spin-1 /2 Kondo 
model known^ to exhibit a KT transition. 

Phase diagrams obtained via NRG are shown in fig. [5] The 
top panel shows the effect of varying the inter-level interac- 
tion U' for fixed U and J^, with behavior that parallels expec- 
tations from the isolated dot limit (sec. Ill Bl and fig. [2}. The 
bottom panel by contrast shows the effect of varying the ex- 
change coupling 7h for fixed U,U', including /h = and an 
antiferromagnetic (AF) /h = —0.5. Note that even for weakly 
AF exchange the USC spin-1 phase still persists, as consid- 
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ered further in sec. HUB I reflecting a ferromagnetic effective 
(RKKY) spin-spin interaction induced on coupling to the lead. 



A. First order transitions on y = x line 

We turn now to the first order transitions permitted by sym- 
metry (sees. (Ill AM Bb ) on the line y = x (£2 — £\ = £ = £ /T). 

To illustrate this, fig.[6]shows the ^-dependence of S; mp (r) 
and rji mp (r), again for fixed U = 20, U' = 7 and /h = 2 (cf 
fig. O, as £ is varied and the transition is approached from 
both sides: £ = e c ± 10~" with n = 4,6,8, 10 ( (b)-(e) respec- 
tively), as well as £ = E c ~ —4.7 itself ((f)). Shown for com- 
parison are the cases ((a)) £ = +5 deep in the Fermi liquid 
regime, with £ here sufficiently large that the degenerate levels 
are barely occupied («j mp ~ 0.25, see top inset); and £ = — 17 
at the p-h symmetric point deep in the USC phase ((g)). 

The stable low-temperature FPs remain of course as before, 
viz. the FI FP for the Fermi liquid phase £ > £ c where the 
global ground state is a singlet; and the USC FP for £ < e c , 
with a doublet ground state. Close to the transition however 
- where the energy separation between these states is tend- 
ing to zero (we denote its magnitude by 7i) - the singlet and 
doublet states will appear effectively degenerate for tempera- 
tures T > ; giving rise in consequence to an entropy plateau 
of Si m p = m 3 seen clearly in fig. [6] with a corresponding 
plateau in the magnetic susceptibility of TXimp = 5 (readily 



3.0 



In 16- 



understood as the mean ((S z ) } = 



imp 

(1x0 + 2 



\) for the 



quasidegenerate states). These are signatures of the 'transi- 
tion fixed point' (TFP), characteristic of systems exhibiting a 
level-crossing transition (see e.g. |[39tl ). On further reducing 
T below ~ 7i, the system is seen to cross over from the TFP 
to one or other of the stable FPs (which crossover in effect 
defines T^). Moreover, as the transition is approached, the 
low-energy scale vanishes - linearly in (e — e c ) as shown 
in fig. [6] (bottom inset), reflecting the level crossing charac- 
ter of the QPT And since T* = precisely at the transition, 
the TFP naturally persists down to T = 0^ 9 (where the ground 
state consists of precisely degenerate global singlet and dou- 
blet states), as evident in case (f) of fig.|6] 

The behavior of the (T = 0) excess impurity charge n; mp is 
also shown in fig. [6] (top inset). In contrast to the continuous 
KT transitions, «; mp is seen to change discontinuously as the 
transition is crossed, commensurate with the first order nature 
of the level-crossing transition. 

A partial progenitor of the latter behavior is in fact appar- 
ent in the trivial non-interacting limit, U = = U' = Ju- Tak- 
ing even (e) and odd (o) combinations of the dot levels 1 and 
2, viz. d ea = (d\a + <ha)/ V2 and d oa = {d\ a - d 2a )/ V2 (cf 
eqn. (O with T,, = T, eqn. (O), only the e-orbital tunnel cou- 
ples directly to the lead and the non-interacting Hamiltonian 
H° reduces to 



H° =\ (£!+ £ 2 ){h e + h ) + V2V(c f ka d ea + h.c. 

k,CT 



(15) 
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FIG. 6: T-dependence of the entropy S; mp (r) (top panel) and spin 
susceptibility T%i mp (T) (bottom) as the first-order transition on the 
line £2 = E] = e is approached and crossed from both sides of the 
transition; for fixed U = 20, U' = 7 and 7h = 2, the transition here 
occurring at e c = -4.73864856302929. Solid lines refer to the USC 
phase (e < e f ), dashed lines to the Fermi liquid (FL) phase. Shown 
for (a) e = +5 deep in the FL phase; (b-e) e = E c ± 10~" with 
n = 4,6,8, 10 respectively; (f) e = e c ; as well as for (g) e = — 17 
at the p-h symmetric point deep in the USC phase. The transition FP 
has a characteristic S; mp = In 3 and TXimp = 51 as indicated; and per- 
sists down to T = precisely at the transition ((f)). Inset, top panel: 
T = impurity charge n imp vs it\ = e. It changes discontinuously 
as the transition is crossed (dashed vertical line), from ;7; mp ~ 0.95 
to 1.31. Inset, bottom panel: The low-energy scale (see text) van- 
ishes linearly in (e — £ c ) as the transition is approached, symptomatic 
of the level-crossing nature of the transition. 



(with i^L the lead Hamiltonian eqn. ([3])). In general, the e and 
o orbitals are coupled by the penultimate term in eqn. ( fT31 ). 
But for the case £2 = £1 = £ of present interest the Hamilto- 
nian is separable, H° = H®+H®, with H® = En a free orbital 
with energy £. The transition in this case thus occurs triv- 
ially for £ = (the p-h symmetric point in the non-interacting 
limit) as the o-orbital - which is unoccupied for £ > - moves 
across the Fermi level, becoming singly occupied precisely 
at £ = and doubly occupied for all £ < 0; such that «i mp 
changes discontinuously from 1 to 3 as £ = is crossed. 

With interactions present the situation is of course much 
less simple. For although the o-orbital remains uncoupled 
from the lead when £2 = £1 , it is then coupled to the e-orbital 
via the non-trivial interaction terms in Hq (eqn. ([T|i), which ac- 
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FIG. 7: Evolution of phase boundaries in the (x, y)-plane for antifer- 
romagnetic J~n < 0. Shown for fixed U = 10 and U' = 5.25 (cf fig.[5] 
bottom), with: J H = -0.55 (long dash line), -0.558 (solid), -0.56 
(short dash) and —0.6 (dotted). Interior regions in each case are the 
USC phase, exterior regions the Fermi liquid. First-order transitions 
are indicated by dots, and occur on the lines y = x and y = —x. 



quire a rather complex (and physically unenlightening) form 
when expressed in terms of e and o operators. We will return 
again to this case in sec. IV B| from the perspective of dynam- 
ics and single-particle 'renormalized levels' . 



B. Weakly antiferromagnetic J H 

As noted above (fig. ID, the USC spin-1 phase survives even 
for weak antiferromagnetic (AF) coupling J\\ < 0, reflecting 
an effective ferromagnetic RKKY interaction induced on cou- 
pling to the lead. For large enough AF |/h I however, the situa- 
tion is clearly different. Only the singlet state of the dot in the 
(1,1) sector is relevant here, and on coupling to the leads one 
expects a global singlet ground state with a stable FI FP which 
is continuously connected to that of the normal FL phases: no 
phase transitions then arise, and the USC phase is eliminated. 

So how is the USC spin-1 phase destroyed as the strength 
of the AF coupling is progressively increased? This is il- 
lustrated in fig. [7] showing phase boundaries in the (x,y)- 
plane for fixed U — 10 and U' — 5.25 (as in fig. |5]l for AF 
/ H = -0.55,-0.558,-0.56 and -0.6. For J H = -0.55, the 
phase boundary has the same form as in fig. [5] consisting of 
the USC phase centred on (x,y)=(0,0), separated from the ex- 
terior Fermi liquid phase by a single boundary line of KT tran- 
sitions except on the line y = x where a first order QPT arises. 

On decreasing /h slightly to —0.558 however, the USC 
phase is seen to split into four distinct domains - symmetric 
as expected under both inversion, and reflection about y = x 
- with the p-h symmetric point in particular now being in the 
FL phase. With a further slight decrease to /h = —0.56, the 
two USC domains straddling y ~x are now eliminated, leav- 
ing two USC regions straddling the line y = — x. This behavior 
persists on further decreasing /h, the remaining USC domains 
diminishing in extent until by /h ~ — 1 they too evaporate and 
the USC phase is eliminated entirely. 

Strikingly, as indicated in fig. [7] one also sees that as the 



USC phase fractionates, first order level-crossing transitions 
arise not only along y = x (as expected on general grounds 
and discussed in sec. HII At : but also along the line y = —x. 

To gain some insight into the above, note that the dif- 
ference in energy (under Hq) between the (1,1) singlet and 
triplet states of the isolated dot is \E$ — Et\ = \Jn\- So for 
|/h| = — Jh ;$ <^(r), and at least close enough to p-h sym- 
metry (x,y)=(0,0) (where «i mp = 2), one expects it necessary 
to include both the (1,1) triplet and singlet states in the low- 
energy dot manifold (higher dot states, such as (2,0), lie con- 
siderably higher in energy provided U' is not within ff(J\\) 
of U). An effective low-energy model within this subspace 
may then be constructed via a Schrieffer- Wolff transforma- 
tion^, the appropriate local unity operator being 1 = lj + Is 
with 1 T = I s; \S = l,S z )(S = l,S z \ and l s = |0,0)(0,0| (with 
\S,S Z ) referring to the (1, 1) triplet or singlet dot states). As 
discussed in Appendix A, the resultant effective model is 



H e s = /i §i • §o + h §2 • §o - I §i • §2 + H L 



(16) 



where as usual S\ and §2 are the spin-1/2 operators for levels 
1 and 2, and §0 = Y.a,a' foo^oa'foo' ls t * le s P m density of the 
conduction channel at the dot (with f^ a = c k a tne cre " 

ation operator for the '0' -orbital of the Wilson chai n 30 ' 31 and 
the number of k-states in the lead). Only exchange scatter- 
ing contributions to eqn. (fT6b are shown explicitly, potential 
scattering being omitted for clarity. The effective exchange 
couplings - viz. the /, > coupling spin i= 1 or 2 to the lead, 
and the direct spin exchange / - are naturally functions of x 
and y; expressions for them are given in Appendix A. 

Eqn . ( fTBT ) is a two-spin Kondo model of the form studied 
in 120J1, so its physics is understoo d 19 ' 20 . A QPT, occurring 
at a critical I c , is obviously driven by the direct exchange /: 
for ferromagnetic AI = I — I c > 0, spins 1 and 2 form a spin- 
1 which is underscreened on coupling to the lead, resulting 
in a residual free spin- 1; while for AF AI = I — I c < by 
contrast, the local singlet Fermi liquid phase naturally arises. 
The resultant QPT is in general of KT form, with one pertinent 
exception^ .: if Ji = Ji in eqn. ( fl6t . then the Hamiltonian is 
separable into distinct singlet and triplet sectors for the spin 
S = §1 +§2, specifically 



eff - 



-# L = It (AS - §0 - ^/) It + |/is :h=h (17) 



(on projecting eqn. ( fTSI l with 1 = It + Is, and using §1 • §2 = 
j(S 2 — |) together with Sis = for all components of S). 

The resultant separability of the Hilbert space for J2 = J\ 
means of course that a first order level-crossing transition can 
occur in this case. As shown in Appendix A (eqn. ( IA14I )). 
this is precisely the situation arising for the present problem 
along (and only along) the lines y = x and y = —x; explaining 
thereby the level-crossing transitions found in fig.|7f 



IV. DYNAMICS AND TRANSPORT 

We now consider dynamics and transport, focussing on the 
T = zero-bias conductance and associated phase shift, 8. 
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A 'Friedel-Luttinger sum rule' for 8 is derived, applicable to 
both the FL and USC spin-1 phases (sec. IIV Bb . A gener- 
alization of Luttinger's integral theorem for the USC phase 
is deduced, and its significant implications for the zero-bias 
conductance determined (sec. IIV Cl >. 



A. Single-particle propagators 

We first summarise basic results for the local single-particle 
propagators, embodied in the 2x2 matrix G(co). Its elements 
Gij(co) are the retarded Green functions for the dot levels, 
i,j G {1,2} (as in sec. [II]), related to the corresponding non- 
interacting propagators G (ft)) by the Dyson equation 



[G(a)}- 1 = [GV)]- 1 -S(o)) 



(18) 



where S(o) is the 2x2 self-energy matrix (with elements 
£y(ft>) = Hfj((o) — iXy(ft))). Using equation of motion meth- 
ods!* 4 ^ the elements of [G°(o)] _1 are given by 



G°(co)] = (ft) + -£;)$,• -r y (fl») (19) 



where co + = (B + iO+, andry(ft)) is the hybridization function 
V,Vi 



r y(ffl) = E 



Tfj{a)-Sij(a>) (20) 



such that Y\ 2 {(0) = rn (a))r 22 (©) (for generality we allow 
here for arbitrary V 2 ,V\). For the standard! flat-band conduc- 
tion spectrum/lead considered (sees. (IIIIIIB ). the imaginary 
part of the hybridization function r^ (co) = Tij (eqns. (1518b ) 
for \co\ < D and zero otherwise; and the corresponding real 
part rfAa) = 0) = at the Fermi level (ft) = 0). 

From eqns. d 1 81 1 9b it follows that G(ffl) has precisely the 
same algebraic structure as G°(ffl), but with r i; (a)) replaced 
by Iy(ft)) defined by: 

fy(ft))=r y (ft))+£y(&)) (21) 
Using eqn. ( fT9l the propagators Gy(ft)) thus follow as 

Gn(ffl) = (ft) + -£ 2 -r 22 (ft))) detG(ft)) (22a) 
G 22 (ft)) = (co + -£i-r n (a))) detG(co) (22b) 
Gia(fl)) = f i 2 (ffl) detG(ft)) (=G 2 i(ft))) (22c) 

with the determinant given explicitly by: 

detG(ft)) = 

[(©+ - ei - f „(©)) (a»+ - e 2 - f 22(a))) - f? 2 (o)] _1 

(23) 

These equations enable the propagators and their spectral den- 
sities Dij(co) = — ilmGy(ft)) to be determined; with self- 
energies obtained in practice via a matrix generalization of 
the standard NRG metho d 34 ' 43 , as outlined in Appendix B and 
discussed further in sec.lVl 



It is also convenient at this point to note an expres- 
sion for the (T = 0) excess impurity charge ni mp ; de- 
fined as the difference in charge of the entire system 
with/without the dot, and hence «; mp = 2-^^Im J^dCO 
[LkGkkM + G„(fi>) + G 22 (co) - L k G£ k (a))] in terms of 
the level propagators G ,-,■(©), the propagators Gkk(ft>) for 
the lead k-states, and their counterparts in the absence 
of the dot, G^ k (co) = [co + — £k]~'- Using equation of 
motion methods it is simple to show that Gkk(ft>) = 
G^(ffl) + Gl k {(0)Z lj V,G ij {(a)VjG" kk ((0), i.e. (via eqn. ®) 
Lk[GkkM - G'^(co)} = -ZijidTijia)/ da>)Gij{co), and 
hence: 

A) 



-1) 



Im 



dco J^Gij(co) 

hi 



8ij- 



dT.jjco) 

dco 



(24) 



For the commonly considered caseA of an infinitely wide flat 
band/lead, where dTjj(co) / dco = for all CO, n\ m? reduces 
to «im P = 2 f_ 00 dft)[Z)ii(ft)) + £> 22 (ft>)] = («i +n 2 ) - i.e. to 
the charge on the dot. For a finite lead bandwidth D (as 
considered here) n; mp is in practice very close to the dot 
charge, although does not coincide identically with it. 

We focus now on the T = zero-bias conductance, given 
from eqn. (0]i by 



G,(r = o) 

Go 



2e 2 ^ 



-T 22 )D ee (co = 0) (25) 



and determined by the Fermi level value of the ee-spectrum, 

(Tn + T 22 )D fe {(o) = EyTy^ImGy(ft)) (eqn. ©); anex- 
plicit expression for which can be obtained using eqns. ( 122123b 
and the CO = behavior of the r, ; (o)). For both the normal 
Fermi liquid phase and the USC phase, the imaginary part of 
the self-energy vanishes at the Fermi level, 



0) = 0. 



(26) 



For the normal FL phase this is of course wholly familiar 1 . 
For the USC phase, we have established it by direct NRG 
calculation of the Ey(©) (it is also consistent with purely 
elastic scattering at the Fermi level for a singular Fermi liq- 
uid^). Given eqn. @, fy(0) follows from eqns. ( 12012 lb 
as fy(0) = -iTy+£ y (0). Using this in eqns. (122123b . and 
defining renormalized single-particle levels in the usual way 1 
by 

E* = £, + l£ (0), (27) 
a simple if tedious calculation gives 

rc(rn+r 22 )Lw(o) = - 



i 

or equivalently 

7r(rn+r 22 )Lw(o) 

with S given explicitly by: 



n 2 



-2r 12 i* (o) 



sin 8 



arctan 



e i*r 22 ■ 



e 2 *r 11 -2r 12 if 2 (o) 



(28) 



(29) 
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B. Friedel-Luttinger sum rule 



The quantity 8 appearing in eqns. ( 128129b is simply the 
static scattering phase shift, given alternatively by 



8 = Imln(detG(0)) 



Imln(detG(fi))) 



I £0=0 



(30) 



(the equivalence of eqns. ( 1291301 ) is readily shown using eqn. 
d23l ): the right hand side of eqn. (f30b also uses arg[detG(© = 
— °°)] = 0, as follows from eqn. d23l together with the fact that 
as | ct> | — > °°, £ 7 (ft)) vanishes while the Ey(co) tend to constants 
(the Hartree contributions). 

We now point out a general result for the phase shift 8. 
From eqns. (1231221 . it follows that 



? ln(detG(0))) = -£^-Afy(a)) 



Integrating this from ffl = — °° to (and noting eqn. (l2Tb ) then 
gives directly from eqn. ( f30b that 



2 "imp 



4 



(31) 



where «; mp is given by eqn. d24l l. and (with Tr denoting a trace) 



Im Tr 



cico — -^^G(co) 



(32) 



is the Luttinger integra l 36 ' 44 (which is dimensionless). 

We emphasise that eqn. OTb is entirely general: applica- 
ble to both the normal Fermi liquid and the USC phases (in- 
deed its derivation does not even require a knowledge of eqn. 
(l26li). For the particular case of the FL phase, Luttinger's the- 



orem gives Ii - 



36,44 



; Ii vanishing order by order in perturba- 



tion theory about the non-interacting limit, reflecting adiabatic 
continuity to the non-interacting limit. In this case eqn. (l3Tb 
reduces to the Friedel sum rulei 3 -^, 5 = fftjmp, relating the 
scattering phase shift to the excess ('displaced') charged in- 
duced on addition of the dot/impurity to the system (and with 
8 £ [0,2ft] for a 2-level dot, since n; mp £ [0,4]). More gener- 
ally, however, 8 and «i mp are related by eqn. (l3"TT l, which we 
refer to as a Friedel-Luttinger sum rule. 

The Luttinger integral for the normal Fermi liquid phase is 
an intrinsic characteristic of it; Il — holding independently 
of the underlying bare model parameters, provided only the 
system is a F L 36 ' 44 . As such, the Luttinger integral is the hall- 
mark of the phase in a rather deep sense. 

The USC spin-1 phase by contrast is a singular Fermi liq- 
uid^. There is no reason here to suppose Ii = 0; and indeed 
it can be shown that the USC phase is not perturbatively con- 
nected to the non-interacting limit of the model. But the obvi- 
ous question arises: as for the FL, does an analogous situation 
arise for the USC phase whereby the Luttinger integral has a 
characteristic value for that phase? 

We answer this question affirmatively, by direct numeri- 
cal calculation (and in several distinct ways). Since the self- 
energies X(co) and Green functions G((o) are calculable from 
NRG, we can calculate Ii directly (eqn. (l32l ) as an co-integral. 



Alternatively, «j mp may be obtained from thermodynamic cal- 
culation (as in sec. [HI) and 8 from calculation of the ee- 
spectrum at the Fermi level alone (as in eqn. d28l ). or alter- 
natively eqn. d29l)); their difference then giving the Luttinger 
integral, //, = 8 — ^«; mp from eqn. (13Tb . We have confirmed 
that the same answer emerges in either way (and for the FL 
phase that II = thereby results). Namely, for any region 
of the (jc,y)-plane where the system is in the USC phase, the 
magnitude of the Luttinger integral is a constant, specifically: 



, , ft 

141 = ^ 



USC 



(33) 



We have repeated the calculations for many different values of 
the bare interaction parameters U, U' and /h. The same result 
emerges; and while the numerics obviously cannot amount to 
a proof, we are confident in the validity of eqn. d33l ). 

Although the magnitude of Ii is constant throughout the 
USC phase, its sign is not. This is a natural consequence of 
symmetry. By considering the symmetries of the propaga- 
tors G(co) = G(co;x,y) and self-energies H(co;x,y) under a 
particle-hole transformation (sec. Ill Ab . it can be shown that 
the Luttinger integral Ii = I^x^y) is odd under inversion, 



h{x,y) = -h{-x,-y). 



(34) 



In addition, as appropriate to the case V2 — V1, the symmetries 
of G, X under the 1-2 transformation (eqn. ( fT3l ) lead to the 
rather obvious invariance under reflection about the line y = x, 

I L (x,y) = I L (y,x) . (35) 

With = f, eqn. (l34l implies the existence of at least one 
bounding curve, of form yt> = f(x) with f(x) = — f{— x), 
across which Ii changes sign from +5 to — 5 ; while eqn. d35b 
for the case V2 = V\ implies that bounding curve to be the line 
y = —x. In practice (by direct calculation, as above) only one 
bounding line is found; and for the case V2 = V\ in particular 
we find: 



h{x,y) 



ft 

2 



y > —x 



= 2 :y< ~ X 



(36a) 
(36b) 



C. Zero-bias conductance 

As above, \Il\ = f is ubiquitous throughout the USC phase, 
as =0 is throughout the normal FL phase. This has im- 
mediate consequences for the behavior of the T = zero-bias 
conductance, given from eqns. (125128131b by 



G C (T = 0) = 2^2 
h 



sin 2 I —n m 



Go h ~ V2 imp 

Since sin 2 (^[«; mp ± 1]) = cos 2 (^«j mp ), it follows that 



G c (T=0) . 7 /ft 

71 = Sln " W'Hmp 

^G V2 



' I 2 "imp 



FL 

USC 



(37) 

(38a) 
(38b) 
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for the FL and USC phases respectively. But as found in sec. 
[Ill] (see e.g. fig. [3j, «i m p varies continuously on crossing the 
line of Kosterlitz-Thouless transitions from the FL to the USC 
phase. Hence from eqn. d38l l it follows that the zero-bias con- 
ductance must jump discontinuously on crossing the QPT; the 
discontinuity on crossing from the FL to the USC phase be- 
ing cos(7T«i mp ), with its sign determined by the value of «i mp 
at the QPT. From direct calculation of single-particle spectra 
we will verify explicitly in sec. [V](see e.g. fig. ITTb that eqns. 
( I38a|38bt are satsified throughout the two phases. Eqn. ( 1381 
is of course equally applicable to the first order level-crossing 
transitions (sec. lIH At . although in this case n; mp itself changes 
discontinuously as the transition is crossed (fig.[6]top inset). 



1.00 



0.75 



<D 

Q ° 0.50 

u 

CNJ 

0.25 



0.00 



1.0 
















0.5 










0.0 










-20 


3 20 








- - , A/ 





-25 -20 -15 -10 -5 5 10 15 20 

co/r 



V. SINGLE-PARTICLE DYNAMICS 

We turn now to ©-dependent single-particle dynamics, here 
focussing primarily on the D ee (a>) spectrum at T = 0. The 
self-energies S(ffl) are obtained from NRG via a generaliza- 
tion of the basic metho d 34 ' 43 to the case of multilevel impuri- 
ties/dots, as outlined in Appendix B. S is thereby calculated 
from 

E(a») = F(a>) [G(co)}- 1 (39) 
where the 2x2 matrix F{(0) has elements 

F^a)) = ((K CT ,4];<» (4°) 

(using conventional notation 1 -^ for the o-dependent 
Fourier transform of a generic retarded correlation function 
((A(/ 1 );J(f 2 )))=-i0(f 1 -f 2 )({A(? 1 ), J §(/ 2 )})); and where 
Hi denotes the interacting part of the dot Hamiltonian, 
given explicitly for the present problem by (eqn. (HJ) 
Hi = l/£ i /l f Jl i | + U'h 1 n 2 — Jn §i -§ 2 . Using the self-energies, 
the fully interacting propagators are then obtained from the 
Dyson equation eqn. (1181) . As for single-level problems, 
calculation of G(co) in this way is numerically stable and 
accurate, and guarantees satisfaction of spectral sum rules, 
J^' oc ,d<X)Dii(co) = l^M 3 - (interleaved NRG/'z-averaging™ 
also being used for optimal calculational accuracy). 

As in sec. [TTT1 we consider explicitly V 2 = Vi', for which 
Fjj = r (eqn. (O), with the ee-spectrum thus given (eqn. (O) 
by: 

D ee {co) = i[Dn(©)+D 22 (ffl)+2Di 2 (©)] (41) 

Fig.[8]shows an 'all scales' overview of 2nrD ee (co) vs co = 
CO/T, for fixed U = 20,U'=7 and J H = 2 (as in fig.[3]for ther- 
modynamics), taking a vertical cut through the (x,y)-phase di- 
agram: the energy of level- 1 is fixed at £j = — \0 — U' = — 17 
{i.e. x = 0), and £ 2 is progressively decreased through the 
Fermi liquid phase towards the transition (e 2c = —6.536...). 

The most important spectral feature is of course the clear 
Kondo resonance straddling the Fermi level. We consider it 
below, but first comment on the qualitative origin of the high- 
energy spectral features, evident most clearly in the three ar- 
rowed peaks shown in fig.[8]for £ 2 = +1. The corresponding 



FIG. 8: Single-particle spectrum 2nTD ee ((a) vs co/T in the Fermi 
liquid phase, for fixed level energy E] = —jU — U' (i.e. x = 0) on 
progressively decreasing £2 towards the QPT (occurring at hi c = 
-6.536..). For U = 20, U' = 7 (i.e. gj = -17) and / H = 2, as in 
fig- EJ with £2 = +1 (solid line), —3 (long dash), —6 (short dash). 
Note that all three cases contain a narrow Kondo resonance strad- 
dling the Fermi level co — 0. For vertical arrows, see text. Inset: 
spectrum for £o = —6.6 on just entering the USC phase, showing the 
absence of a Kondo resonance. 

evolution of «; mp vs e 2 is shown in fig. [3] (inset), from which 
it is seen than «j mp ~ 1.2 for e 2 = +1 - sufficiently close to 
unity that we can interpret the high energy spectral features 
as removal or addition excitations from the singly-occupied 
(«i,« 2 )=(l,0) state of the isolated dot. The removal exci- 
tation from dot level-1, contributing as such to the D\i((o) 
constituent of D ee ((o) (eqn. (RTT)). thus correponds trivially to 
£z)(l,0) — £0(0,0) = £1 (in the notation of sec. Ill Bt . i.e. lies 
below the Fermi level at (co/T =) &> = ej = — 17 here, gener- 
ating the lower 'Hubbard satellite' seen clearly in fig. [8) its po- 
sition, dependent at this crude level of description only on t\ , 
varies only slightly on further decreasing e 2 in the FL regime. 

Two addition excitations lying above the Fermi level are 
also seen in fig. [8] for £ 2 = +1 (arrowed). The lowest cor- 
responds to electron addition to level-1, and hence shows up 
(again via fln(fi))) as an excitation at £d(2,0) — £0(1,0) = 
£1 + U; thus lying at <» = <a + = E\ + U = +3 as seen in the 
figure. The second excitation corresponds to addition to level- 
2, contributes as such to the D 22 (fi)) constituent of D ee (co), 
and thus corresponds to -Ed(1, 1) — -Ed (1,0). Since there are 
two distinct (1,1) dot states - triplet and singlet - two such 
excitations in principle arise; separated in energy by /h and 
occurring at 6) = G>t = £2 + U' — ^Jh (for triplet (1,1)) and 
co = cb s = e 2 + U' + |/h for the singlet. As evident from the 
figure, coupling to the leads in practice blurs these excitations, 
so that only a single spectral feature is seen. 

On decreasing e 2 from +1 the 0>s/t excitations (which as 
above depend on e 2 ) decrease, and become comparable in en- 
ergy to the co + (= +3) excitation; so that as seen in fig.[8]the 
high-energy addition excitations in practice merge to a form 
a single peak, which on decreasing e 2 through the FL phase 
moves towards - but does not reach - the Fermi level. In- 
stead, the single-particle spectrum in the immediate vicinity 
of the Fermi level co = is naturally dominated by the narrow 
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FIG. 9: For same parameters as fig. [8] showing a close-up of the 
Kondo resonance on approaching the transition (at £o c = —6.536..) 
from the FL side: 27tFD ee (co) vs a>/Y for &> = —6.1 (solid line) and 
—6.2 (long dash). The Kondo resonance collapses 'on the spot' as the 
QPT is approached and the Kondo scale Tk — > 0. The short dashed 
line shows the spectrum for ft = —6.54 on just entering the USC 
phase; it is featureless on these scales, with no Kondo resonance. 
Inset: Scaling of the Kondo resonance on approaching the QPT from 
the FL side. Both FL spectra in the main figure collapse to a universal 
scaling resonance as a function of Co/Tk (their individual T%s differ 
by more than an order of magnitude). 



low-energy Kondo resonance evident in fig. [8] 

The evolution of the Kondo resonance itself is shown in 
close-up in fig. [9] As £2 — ► £2c+ from the FL side, it narrows 
progressively - reflecting the incipient vanishing of the Kondo 
scale Tk known from thermodynamics (sec. [Till figs. d3l41 >) - 
and collapses 'on the spot' at the transition itself, where Tk 
vanishes. As a corollary, in the USC phase just on the other 
side of the transition the Kondo resonance is simply absent; 
as seen in fig. [9] (for £2 = —6.54) where the USC spectrum 
is constant on the low &> = (o/Y scales shown. The inset to 
fig. [8] also shows this USC spectrum on an 'all scales' level, 
showing that while the Kondo resonance is absent here, the 
high-energy features discussed above evolve in a smooth way 
from those arising in the FL phase. 

Since the Kondo scale Tk vanishes as the QPT is ap- 
proached from the FL side, one expects the Kondo resonance 
to exhibit universal scaling in terms of it. That this is so is seen 
in fig.|9](inset), where both FL spectra shown in the main fig- 
ure collapse to a universal scaling form as a function of Co/Tk- 
Note also that while we have scaled the spectra here in terms 
of Tk obtained from S; mp (r) (as in sec. UTIt. we could equally 
have defined Tk spectrally - e.g. via the width of the Kondo 
resonance - and likewise obtained universal scaling behavior. 
The essential point is simply that there is only one vanish- 
ing low-energy scale as the QPT is approached, and different 
practical definitions of it are all fundamentally equivalent. 

The subsequent evolution of the spectrum in the USC phase 
is shown in fig.[T0] Not far into the USC phase (£2 = —7), the 
spectrum lacks a Kondo resonance, as above. However on fur- 
ther decreasing £2, a second Kondo resonance straddling the 
Fermi level is seen to arise. It is in fact well developed al- 
ready by £2 = —10, and narrows progressively as £2 decreases 



Q 

U 

CM 



1.00 



0.75 



0.50 



0.25 



0.00 




FIG. 10: USC phase. For same parameters as figs. J8l9b . 27tFD ee (co) 
vs co/Y for £2 = —7 (short dash line), —10 (dashed) and the center 
of the USC phase £2 = — 17 (p-h symmetric point, solid line). The 
Kondo resonance which develops here is that for a spin-1 Kondo 
model — ; see also sec. IV Al (fig.ll3b. 



towards the center of the USC phase at the /7-/j-symmetric 
point £2 = — 17. The origin of this behavior is readily guessed 
from «i mp vs £2 (fig. [3] inset). For although the transition it- 
self corresponds to a 'mixed valent' «; mp w 1.4, on entering 
the USC phase «; mp increases quite rapidly; such that even 
by £2 = —10, «; mp is close to 2. Here one expects the sys- 
tem at low energies to be described asymptotically by a spin-1 
Kondo model, and hence the second Kondo resonance to be of 
that ilk. This is indeed so; we discuss it further in the context 
of fig. [13] below. High-energy spectral features in this regime 
are also naturally interpretable in terms of single-electron ex- 
citations to/from the (1,1) triplet ground state of the isolated 
dot; e.g. at the p-h symmetric point, all addition/removal ex- 
citations to/from both levels 1 and 2 have the same magni- 
tude, \e\ + U' — \fn\, giving rise to the symmetrically dis- 
posed Hubbard satellites at |<5| ~ 10.5 seen in fig. ITOl 

Finally and importantly, fig. QT| verifies the predictions of 
sec. IIVCI for the behavior of the zero-bias conductance on 
crossing the QPT. The Fermi level spectrum 2nYD ee (Q) vs £2 
is shown in both phases, and compared explicitly to eqn. (1381 
with «; mp obtained from an independent thermodynamic NRG 
calculation; the agreement being excellent. 



A. Kondo antiresonances 

In the example considered above the QPT is associated with 
a collapsing Kondo resonance in the FL phase; and hence nat- 
urally with a decrease in the zero-bias conductance on cross- 
ing into the USC phase. From eqn. d38l ), the latter behav- 
ior is generic provided ni mp at the transition lies in the in 
terval n 



£ [5,5] (by symmetry eqn. ( fT4b we can consider 



"imp € [0,2] rather than the full range [0,4]). If however n^n at 



the QPT lies in the range [|,2] or [0, 5], then eqn. d38l predicts 
generically an increase in the conductance on crossing from 
the FL to the USC phase. One might thus intuitively expect 
such behavior to be associated with a vanishing Kondo an- 
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FIG. 11: 27tTD ee (co = 0) (equivalently the zero-bias conductance, 
see eqns. | |25I38| |) £2 in both phases on either side of the Kostelitz- 
Thouless transition (dashed vertical line at &i c = —6.536.., same bare 
parameters as figs. f8l-U0b). Crosses show the co = spectra deter- 
mined from NRG, while solid lines show sin 2 (-| "imp) (in FL phase) 
and cos 2 (^H; mp ) (in USC phase) with nj mp obtained from a thermo- 
dynamic NRG calculation; verifying the predictions of eqn. d 3 8 1 > - At 
the QPT, «[ mp = 1.43, whence the spectrum/conductance decreases 
discontinuously on crossing from the FL to the USC phase, «; mp it- 
self evolving continuously (fig.[3]inset). 



indeed contains a Kondo antiresonance in the FL phase, here 
with 2^rD ?( . {co = 0) = throughout. This antiresonance like- 
wise vanishes on the spot as the transition is approached and 
the Kondo scale Tk — ► 0; and as it does so exhibits scaling 
as a function of Co/Tk (fig. [T2] right inset), the low-frequency 
spectral behavior being 2nTD ee {co) °= (co/Tr) 2 , symptomatic 
of a normal Fermi liquid. 

Note that the general predictions of sec.||V]are neatly exem- 
plified by the above results: since «; mp = 2 everywhere along 
the y = -x line, eqns. ( 125138b yield 2nTD ee {0) = 1 in the USC 
phase and in the FL phase (as confirmed by direct calcula- 
tion, figs. ( I12I13H : and hence that the zero-bias conductance 
G C {T =Q)/Gq in this case increases by precisely the conduc- 
tance quantum 2e 2 /h on crossing the QPT into the USC phase. 

Fig. [T3l continues fig. [12] into the USC phase, showing the 
ee-spectra for y = 5,3 and 0. As for its counterpart in fig. 
[TUl the Kondo resonance which develops in the USC phase 
is that for a spin-1 Kondo model 2 ^. As shown in Appendix 
A, its low-energy scale Tg varies with the bare interaction 
parameters as (modulo an immaterial prefactor) 



exp 



(U- 



v 2 



U 



(42) 



1.50 
1.25 
1.00 



1.0 
0.5 
0.0 







1.0 








0.5 


\ /' 
\ / 






A 0.0 


. v . 




FIG. 12: 27zTD ee {co) vs co/T along the line y = — x (for same in- 
teraction parameters as figs. ([81-11 It), decreasing y = £2 + \ U + U' 
through the QPT (at y c = 6.36..) from the FL side y > y c . Shown for 
y = 6.55 (short dash) and 6.40 (long dash) in the FL phase, and y = 
6.30 (solid) just into the USC phase. A clear Kondo antiresonance 
at the Fermi level in the FL phase is seen (with 27tTD ee {0) = 0). It 
collapses 'on the spot' as Tk —* and the QPT is approached; and 
exhibits scaling as a function of Co/Tk as it does so (right inset). Left 
inset: as main figure, on an expanded co/T scale. 



Hence on decreasing y through the USC phase, the Kondo 
resonance becomes increasingly narrow as Tg decreases to- 
wards its smallest (but non-zero) value occurring at the p-h 
symmetric point y = (= x); and as shown in fig.[T3] univer- 
sal spectral scaling as a function of co/Tjt =l thereby arises. 

Fig.[l3](inset) also shows the clear cusp-like behavior of the 
spin- 1 Kondo resonance as — > 0, known from study of the 
spin-1 Kondo model itself (with spectra inferred from the t- 
matrix of the Kondo model). This behavior is characteristic of 
the singular Fermi liquid 2 - nature of the underscreened spin- 
1 phase; specifically the weak ferromagnetic coupling of the 
residual spin- j to the metallic lead, resulting in logarithmic 
corrections to Fermi liquid behavior. As |o|/7^ =1 — > we 
find 



2nTD ee {co = 0) - 1 



ln 2 (|to|/^ 



S=l N 



(43) 



(with b a constant), the leading logarithmic correction here 
stemming from the leading low-O) behavior of the self- 
energies LUw) ~ l/ln 2 !!©!/^ 1 ); and which form is in 
agreement with that of 112111 for the spin- 1 Kondo model itself. 



tiresonance as the transition is approached from the FL side. 

That this indeed arises^ is illustrated in fig. [12] where dy- 
namics on the line y = — x are considered; along which, by 
symmetry (eqn. (fPfl) ). «j mp = 2 regardless of phase (the spec- 
tra are likewise readily shown to be symmetric in co). For bare 
interaction parameters U = 20, U' = 7 and /h = 2, we de- 
crease (y/r =) y = £2 + jU +U' (eqn. (fTTT >> across the transi- 
tion occuring at the critical y c = 6.36.., from the FL side (y > 
y c ) to the USC phase. As shown in the main figure D ee (co) 



B. The y = x line 

As considered in sec. IIII Al in regard to thermodynamics, 
the transition occurring along the line £2 = £\ {i.e. y = x) is a 
first order level-crossing transition, as permitted by symmetry 
for V2 = V\ . Here we consider the £2 = £1 line again, from the 
perspective of dynamics, and the resultant channel separabil- 
ity arising in the 'even/odd' representation as now discussed. 
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FIG. 13: Continuing fie. 1 121 into the USC phase along the line y = 
—x: 27lTD ee (co) vs Co/T for y = 5 (short dash), 3 (dash) and the p-h 
symmetric point y = = x (solid). The Kondo resonance developing 
in the USC phase is that for the spin- 1 Kondo model!!.; and (inset) the 
three spectra in the main figure show low-energy universal scaling as 



a function of ft>/7£ 



with t£~ 



the spin-1 Kondo scale (see text). 



1. Even/odd basis 



In previous sections the elements of the Green function ma- 
trix G(co) have been considered as the propagators for the dot 
levels, viz. Gy(ft)) with i,j G {1,2}. Equally, one can take 
even/odd combinations of the dot levels, viz* deo — 
d2a)/V2 and d oa = [d\ c — ctf 2 <T)/v2, and consider G((o) in 
an e/o representation; with elements G a p (ft)) given explicitly 

by Gee((o) = \ [G n (ft)) +G 22 (ft)) ±2Gi 2 (ft))], with G eo (ft)) = 

j[Gu(o>) — G 22 (ft))] (= G oe (ft))) for the off-diagonal ele- 
ments. For £ 2 7^ £i in general, there is no particular advan- 
tage in working with the e/o representation. However along 
the line £ 2 = £i where levels 1 and 2 are equivalent by symme- 
try, Gn (ft)) = G 22 (ft)) and hence the off-diagonal G eo (co) = 0. 
G(co) in the e/o representation is then purely diagonal for all 
ft), with elements: 



Gee(G>) = Gn (ft)) ± Gi 2 (ft)) 

Using eqns. ( 12212312 lT ). one obtains 

G ee (co) = [ft)+-e-2r(ft))-£ ee (ft))] 



G o0 (ft)) 



co 1 



(44) 

(45a) 
(45b) 



where £ = £i = £ 2 denotes the common level energy, the hy- 
bridization function is r(co) (= r i; (ft)), eqn. ( f20l > with V 2 = 
Vi), and the ee/oo self-energies are given simply by: 



Zee((0) = En(ft)) ± £i 2 (ft)) 



(46) 



Notice from eqn. (I45bl i that there is no direct hybridization 
(r(ft))) contribution to G ou (co), reflecting the fact (sec. IIII Ab 
that for £ 2 = £j the o-orbital is not directly coupled to the 
lead. In the non-interacting limit the o-level is thus entirely 
free, G^,(ft)) = [ft) + — e] , but in general the o/e levels are 
coupled via interactions, as embodied in £ „(ft)) ^ 0. 



Since G(a>) is diagonal in the e/o representation, 
detG(ft)) = G ee ((0)G oo ((0), and hence from eqn. ( f30b the 
static phase shift 8 is separable into e and o channels, 

8 = 8 e + S . (47) 

A short calculation using eqn. d30l ) (together with YJ aa (co = 
0) = from eqn. d26l l) then gives 



Se 

8 



arctan 



2r 

OH 



arc 1311 ! ) = 7ld (- £ o) 



(48a) 
(48b) 



where each 8 a <5 [0, n], Y (= r (ft) = 0)) is the usual hy- 
bridization strength (eqn. (0), and 9(u) is the unit step func- 
tion. The £^ denote the renormalized e/o levels, given by (c/ 
eqn. d27| )) 



■4«(o) 



(49) 



with £* = £^(x) such that £^(x) = -£* (— x) (via a trans- 
formation, sec. Ill Al l. Likewise, considering <5lnG aa (ft))/(9ft) 
and repeating the calculation leading to eqn. OTT ). gives 



8a — 'J "imp, a 



where (cf eqn. ( l32l ) 

/« = Im / 

J — c 



<9E aa (ft)) 
c/co — 5 G aa (co) 

(7ft) 



(50) 



(51) 



is a Luttinger integral for channel a = e or o (with = /" (x) 
such that (r) = —I" ( — x) under inversion); and «i mp , a is the 
excess impurity charge associated with channel a, given by 

dT(co 



nmp.e 



2i^Im / d(OG ee ((0) 



1-2- 



(9co 



»im P ,o = 2 ^Im / do) G 00 (ft) 



(52a) 
(52b) 



such that the overall n 



t imp,o: 



imp " 



I imp,e 



nmp,( 



„ (eqn. d24b). and with 
(— x). Since the behavior of relevant 
x is as specified above, we 
can focus oni = £ + jU + U' > 0; and do so in the following. 



'imp. a 



quantities under inversion x 



2. Results 

The charges «i m p,a may be calculated directly from NRG. 
Their evolution with x = x/Y is illustrated in fig. [T4l (inset), 
on decreasing e = e/Y with U,U', /h fixed at the same values 
used in fig. [6] for thermodynamics; the level-crossing transi- 
tion here occurring at x c = 12.26.. (£ c = —4.73..). On cross- 
ing the transition from the FL side (x > x c ) to the USC phase, 
«i mp increases discontinuously from - found to be its value 
for all x > x c — to «imp,o = 1, which constant value is like- 
wise found throughout the USC phase x < x c . For the e- 
channel by contrast «i mp e is not fixed in either phase; but it too 
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FIG. 14: Renormalized odd level e* /T (eqns. <49l54t ) vs x = x/Y = 
£ + jU + U' (for fixed = 20, U' = 7 and J H = 2). The level- 
crossing transition at x c = 12.26.. is marked by an arrow; e*, ~ 
x — x c = £ — £ c vanishes linearly as the transition is approached. 



Inset: Behavior of «i mp , (solid line) and n u 



(dashed) on cross- 



ing from the FL phase (x > x c ) to the USC phase. Both jump dis- 
continuously on crossing the transition FL — * USC, «i m p !0 upwards 
from to 1 and n[ mp e downwards as indicated; such that the overall 



eqn. ( IB5b in the diagonal e/o representation]. The generic 
i-dependence of e* = e* /T is illustrated in fig. [T4j It evolves 
continuously for all x > (the divergence on approaching the 
p-h symmetric point x = at the center of the USC phase re- 
flects via eqn. d54]i the fact that F^(0) — ► — 1 T as x — > 0*). 
In particular, in the USC phase < x < x c , the renormalized 
level e* < 0, while fori > x c in the FL phase, e* > 0; the level 
vanishing linearly as the QPT is crossed, 

e ~ x — x c = e — e c . (55) 

And for large enough £> 1, where both the e and o levels are 
in practice empty and interaction effects embodied in £ oc are 
thus irrelevant, e* — > e (the 'bare' level energy, see eqn. d49l)). 

The above results then enable the o-channel Luttinger inte- 
gral I'l (eqn. (f5T)) to be deduced. Since £* < [> 0] in the 
USC [FL] phase, eqn. ( I48bl i gives a phase shift 8 () = n in the 
USC phase < x < x c , and 8 () = in the FL phase x > x c . 
Combining this with eqn. d53l ) for n imv o , the Luttinger inte- 
gral I'[ = 8„ — f rtjmp )0 (eqn. d50l l) follows directly as 



^lmp 



o + ' ! imp e increases as shown in fig.|6](top inset). 





K 

2 



: FL, x > x c (56a) 
: USC, 0<x<x c (56b) 



jumps discontinuously, decreasing as the FL— ^USC transition 
is crossed. The latter behavior is physically natural, since 
the piling of charge into the o-orbital which accompanies the 
transition increases Coulomb repulsions with electrons in the 
e-orbital, which the concomitant reduction in rii mp e acts to 
offset. The two effects do not however cancel, the overall 



- which result we have also verified by direct computation of 
I[ itself, eqn. dBTT l. 

For the e-channel by contrast, direct calculation of I[ gives 
I[ = in both the FL phase and the USC phase, 



II = 



FL and USC 



(57) 



nmp 



nmp.o 



«im P , e (shown in fig. |6] top inset) increasing as The total Luttinger integral I L = l° L 



the transition is crossed. The behavior just described is redo- 
lent of, but distinct from, that occurring in the non-interacting 
limit discussed in sec. IHI Al where at the transition, in that 
case occuring for e c = 0, «i mp , jumps discontinuously from 
to 2 but with no concomitant change in «i mp e since there are 
no interactions present. Moreover since the transition is gener- 
ically accompanied by occupancy of the o-orbital, one intu- 
itively expects the requisite critical e c for the transition with 
interactions present to be reduced below its non-interacting 
counterpart e c = 0, in order to offset the increased interac- 
tions; as indeed is found. We also add that the behavior found 
is not specfic to the interaction parameters used for illustra- 
tion; in particular that 



I[ thus vanishes as re- 



quire d 36 i 44 throughout the FL phase; while for the USC phase 
eqns. ( I56bl57l ) agree as they must with the general result eqn. 
d36at for Ii (x,y) (which is not confined to the v = x line). Note 
further, using eqn. d57| i, that eqns. d48a|50l ) give 



e* = 2T tan (f [l - « imp , ? ] ) : FL and USC 



(58) 



independently of the phase (as again verified by separate cal- 
culation of £* and «i m p,f). From the x-dependence of «i mPjf . 



hmp,o 





1 



FL, x > x c 
USC, x < x c 



(53a) 
(53b) 



is found to occur ubiquitously. 

We consider now the static renormalized levels, calculable 
from eqn. (149b ; or, for £*, equivalently from 



(54) 



1+F*(o> = 0) 



where F oo {(0) = F n (co) - F 12 (o)) (with F ( f f} (co) its real part), 
and the Fjj(a>) are given by eqn. (|40b and calculated directly 
via NRG [eqn. ( l54l l follows from eqn. d49l together with 



illustrated in fig. [L4] (inset), eqn. (|58l shows that progres- 
sively decreases as x is decreased through the FL phase, in- 
creases discontinuously as the FL^USC transition is crossed, 
and in the USC phase decreases monotonically as x is de- 
creased towards the p-h symmetric point x = 0, where «i mPje = 
1 and hence e* = (and with e* (x) for x < following from 
the symmetry E*(x) = — £*(— x)). 

As a brief illustration of single-particle dynamics along the 
y = x line, fig. [TBI shows the evolution of the o-orbital spec- 
trum D O0 ((0) = — jImG m ,(ffl) on decreasing.? through the FL 
phase (main panel), across the transition and into the USC 
phase (right inset). In the vicinity of the QPT coming from 
the FL side, a strong low-frequency spectral resonance (for 
ft) > 0) is seen to develop; becoming a pole at the Fermi level 
precisely at the transition, and crossing smoothly to ft) < 
in the USC phase. The position of the resonance tracks the 
vanishing renormalized level e* (eqn. d55ll), the ft) = pole 
at the transition reflecting e* = (from eqn. (I45bb . using 
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FIG. 15: o-orbital spectrum TD 00 (co) vs co/T on decreasing x = x/T 
in the FL phase (for same interaction parameters as fig. [14] with the 
critical x c = 12.26..). Shown for x = 18, 17, 16, 15, 14, 13 and 12.5 
(from right to left). Right inset: as main figure, but in the USC phase 
for x = 12 (dotted line), 7 (dashed) and at the p-h symmetric point 
x = (solid). The o-spectrum evolves continuously as the transition 
is crossed, a pole occuring at the Fermi level precisely at the transi- 
tion where e* vanishes. Left inset: e-orbital spectrum 2nTD ee (co) vs 
ffl/rjust on either side of the QPT, forx= 12.261351440 in the FL 
phase (solid line) and x = 12.26135 1420 in the USC phase (dashed). 
Here the entire spectrum changes abruptly on crossing the transition. 



£ 00 (co = 0) = £« o (0) together with eqn. (O, the Fermi level 
spectrum is given generally by D 00 (co = 0) = 5(e*)). In the 
vicinity of the transition, the renormalized level e* is the coun- 
terpart of the low-energy scale introduced in sec. IIII Al in 
respect of thermodynamics (see e.g. fig. |5); and e* both 
vanishing linearly as the transition is approached, and control- 
ling the low-energy behavior of appropriate thermodynamics 
and single-particle dynamics respectively. 

We also add that, as expected on physical grounds, the 
vanishing o-orbital renormalized level does not show up in 
the corresponding e-channel spectrum D ee (co), which as illus- 
trated in fig.[l5](left inset) changes in a wholly discontinuous 
fashion on crossing the transition; commensurate with the in- 
herently first-order nature of the transition along the y — x line. 



VI. EXPERIMENT 

We now consider the experiments of Kogan et al^ on a 
GaAs-based single-electron transistor at low temperature (T), 
embodied in the differential conductance as a function of gate 
voltage AVg (measured relative to a reference voltage), and 
also the bias (or source-drain) voltage V sd . On varying the gate 
voltage, the resultant conductance maps shown e.g. in fig. 1 of 
ifTTll (see also the theoretical fig. [l6]below) show clear zero- 
bias Kondo peaks arising in the centers of adjacent Coulomb 
blockade valleys; one valley thus being associated with an odd 
number of dot electrons and the other with an even number. 
The former valley, which extends over a relatively wide AV g 
range, is naturally interpreted!! as the normal FL, or 'singlet 



phase'; while the latter, extending over a narrower AVg range, 
is interpreted!], as the 'triplet phase' (i.e. the USC phase). 
In considering theoretically the conductance, 



G c (r,v sd = o) 

(2e 2 /h)G 



dco 2nYD ee ((o) (59) 

; dCO 



is exact^Z at zero-bias (as before we consider explicitly V2 = 
Vi), with D ee (co) the spectrum at the temperature of interest. 
At finite bias by contrast, nothing exact can be said with the 
methods at hand. To treat approximately V sd 7^ we neglect 
explicit dependence of the self-energies on V sd . With this stan- 
dard approximation G c = G c (r, V sd ) is readily shown to be 



G c (r,v sd ) 

(2e 2 /h)G Q 



dco 



-df L (co) , -df R (co) 



dco 



dco 



{f {CO) = [e«>l T _ 
rectly to eqn 

G c (7/ = 0,V sd ) 
(2e 2 /h)G 



2nTD ee (co) 
(60) 

where f v {co) = f(co ± \eV s d) for lead v ~R/L respectively 
1]~'). For V sd = 0, eqn. ( f60b reduces cor- 
; while for T = it yields 

tzF [D ee (co = ±eV sd ) +D ee (a = -^eV sd )] 

(61) 

in terms of the T = spectra. In the above we have taken 
a symmetric voltage split between the R/L leads. From eqn. 
( |60t this gives G c (r,V sd ) = G C (T, — V sd ), which symmetry is 
rather well satisfied in experiment (fig. 1, ifTTll ). 

Under application of a gate voltage, the level energy E\ °^ 
AVg, and one expects the level spacing Ae = £2 — £1 to be es- 
sentially fixed^l The experimental 'trajectory' in the (£i,£2) 
(or (x,y)) plane upon varying AVg is then as indicated schemat- 
ically in fig. 12; a), viz. y = x + Ae (i.e. £2 = £1 + Ae). In 
this regard an interesting symmetry arises. Indicating explic- 
itly the x,y dependence of D ee {co) = D ee (co;x,y), it is readily 
shown that under the p-h and 1-2 transformations of sec. Ill Al 
D ee (co;x,y) — D ee (—co;—y,—x). Employing this in eqn. d60b . 
noting that [dfi(co)/ dco + dfn(co)/ dco] is even in CO, gives 



Gc(T,V s(i ;x,y) = G c (T,Vsd;-y,-x) 



(62) 



for the conductance G c (r,V s d) = G C (T, V s d',x,y). That is, 
the conductance is symmetric under reflection about the line 
y = — x. Now the 'trajectory' y = x + Ae is perpendicu- 



-x, and intersects it for x 



\Ae, i.e. 



■hu- 



tf) for £1 =-\Ae-\U -U' = £i, n 



lar to the line y 

(since x = e\ , ^ ' " 1 fc i — 2^ ~ 2 1 - 
Since the phase boundaries are also symmetric under reflec- 
tion about the line y = —x (figs. (I2I5I7H , this value £j = £j m 
- and hence the corresponding AV 2 , m (£1 « AV g ) - is thus 
the midpoint of the triplet phase; such that the conductance 
should be an even function of £1 — £1 m , or equivalently of 
AVg — AVg. m . This symmetry is quite well satisfied in ex- 
periment (and is of course obeyed precisely in the theoreti- 
cal results). Fig. [16] (middle panel) shows the experimental 
zero-bias conductanc e) 1 *i 46 (crosses) as a function of AV 2 (in 
mV), together with corresponding theoretical results (as de- 
tailed below). The midpoint of the triplet (T) phase is readily 
identified as AV £ 



10 mV, about which the experimental 
conductance is indeed seen to be quite symmetric. And the 
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FIG. 16: Middle panel: Experimental zero-bias conductanc e) 11 ' 46 
(crosses) vs AVg (in mV). Theoretical results obtained as described 
in text are also shown, for T /T = (solid line), 5 x 10~ 3 (long dash) 
and 10~ 2 (short dash); for T = 0.5meV these correspond respectively 
to T = 0, 30mK and 60mK. The positions of the theoretical phase 
boundaries (T = 0) between the USC triplet (T) and Fermi liquid 
singlet (S) phases are as indicated. Upper panel: Kondo scales de- 
termined as in sec. |ni] and shown as log(r^-/r) vs AVg (for the T 
phase, T% = T^ 1 ). Bottom panel: Theoretical differential conduc- 
tance map in the (V s[ j, AV g )-plane, with the grey-scale code indicated 
(in units of 2e 2 /h). Shown for T = 30mK, choosing T = 0.5meV 
(see text). The experimental counterpart is shown in fig.l of fiUl . 



experimental conductance map shown in fig. 1 of JTlll (cf fig. 
[TrSlbottom panel) is also clearly rather symmetric about AV" g m . 

To compare directly to experiment we must specify the di- 
mensionless interactions ,0' ,1a (eqn. ([Toll). Ae = Ae/T, the 
relation between £\ = £\ /T and AVg, and finally the hybridiza- 
tion strength F. This is obviously a large parameter space, and 
our intent here is simply to employ what we regard as a rea- 
sonable set of 'bare' parameters. For a typical dot the relative 
hierarchy of energies satisfies^ |/fj| <C Ae -C 0; with which 
the specific parameters we use here concur, |/h| = 0.5,Ae = 
4.5 and = 12 (and with and Ae in excess of unity, con- 
sistent with the occurrence of charge quantization towards the 
centers of the Coulomb blockade valleys). No attempt to ex- 



plain experiment on the assumption 0' = was found to be 
successful, even qualitatively, on varying the bare parameters. 
The main reason (as evident from inspection e.g. of fig. |2c) 
or fig.[5](top, (c)) is that the resultant width of the T phase (in 
£i or AVg) is much too large compared to that of the singlet 
(S) phase, and as such not qualitatively consistent with exper- 
iment^ (fig. [Tol l. For the results shown here we have used 
0' = 6 = 0/2 (although tolerable variations from this value 
give comparable agreement with experiment). 

From the discussion above, the relation between AVg and 
£i is of form AV g = c[E\ — £j m ] + AV~ g . m where the propor- 
tionality constant c is to be determined (as above, AV g . m = 
— 10m V and ei >m = — \Ae — 10 — 0'). For a chosen set of 
,0' ,/h,A£, the theoretical zero-bias conductance at T = 
is calculated from eqn. d59l ) as a function of £\. It is then 
scaled onto the experimental results shown in fig. [16] (mid- 
dle), over the AV g range above ~ 35 — 40m V. We choose this 
range because here the system is beginning the approach to the 
'empty orbital' regime of «j mp <C 1, where one does not expect 
any appreciable T-dependence to the conductance [the exper- 
imental T is not known with certainty 4 ^ for although the ex- 
periments were performed at the refrigerator base temperature 
of ~ 12mK, the electron temperature, T, was not determined; 
although it is believed to be < 40mK 46 ]. With this procedure 
we determine the constant c, which is then fixed and used for 
all AVg (and T); as well as the dimensionless constant Go re- 
flecting (sec. |llj the relative asymmetry in tunnel coupling to 
the leads (from scaling the vertical axis in fig.HS] and leading 
to Go — 0.8 - as is obviously reasonable even from cursory 
inspection of the experimental data). 

In comparing to experiment, an obvious key element is the 
relative widths (in £i or AV g ) of the S and T phases, the for- 
mer being considerably wider than the latter in experiment. 
This we naturally find to be influenced significantly by the ex- 
change /h (and to a lesser extent by Ae), which is optimised 
accordingly. For the results shown here, we find /h = —0.5 
to be optimal. Its magnitude is small, as expected, although 
its sign is antiferromagnetic. This is not however unreason- 
able, for on coupling to the leads as mentioned in sec. [Ill] 
an AF bare /h still generates an effective ferromagnetic spin- 
coupling via an RKKY interaction (as evident in the very ex- 
istence of the USC triplet phase for weakly AF bare /h, and 
which effect is in fact largest for the case Vz = V\ we consider 
explicitly). 

With the above we calculate the T — zero-bias conduc- 
tance, shown in the middle panel (solid line) of fig. [16] for 
U = 12,0' = 6,A£ = 4.5, fa = -0.5, with the two T/S phase 
boundaries indicated in the figure. The T phase, symmetri- 
cally disposed about the midpoint AV g , m = — lOmV, occurs 
in the interval — 13.45 mV < AVg < — 6.55mV (correspond- 
ing to —15.2 < £i < —13.3); with «; mp at the phase bound- 
aries of «; mp = 1 .87 (upper boundary at AV g = — 6.55mV) and 
«imp = 4— 1.87 = 2. 13 (lower boundary), such that in accor- 
dance with eqn. d38l of sec. IIVCI the zero-bias conductance 
increases on crossing from the S (FL) to the T (USC) phase. 

The resultant Kondo scales as a function of AV g are shown 
in the top panel of fig. [16] (obtained as specified in sec. [HI] and 
with Tk = 7^ =1 for the T phase). Since Tk vanishes as the 
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QPT is approached from the S (FL) side, finite-temperature 
effects will obviously be most significant in the vicinity of 
the transition. Fig. Q~6] thus shows the zero-bias conductance 
at two non-zero temperatures, T /T = 0.005 and 0.01. While 
there is not much net difference between the two, each has the 
effect of significantly increasing the conductance in the vicin- 
ity of the transition, and leads to what we regard as rather good 
overall agreement with experiment. Over the r -range shown 
the conductance 'inside' the T phase does not erode as rapidly 
as one might like, reflecting the fact that T% = Tg therein is 
in excess of the 7"s shown; although one could likely improve 
on this with a bare parameter set for which 7^ =1 inside the T 
phase is somewhat smaller. The temperature range considered 
here is also entirely reasonable in relation to the experimen- 
tal T discussed abov o 11 ' 46 ; with T = 0.5meV (as employed 
below) the temperatures shown correspond to T = 30mK and 
60mK respectively. 



A. Conductance maps 

The bare parameters specified above are fixed. Using eqn. 
( |60b the differential conductance, as a function of gate and 
bias voltages, may now be calculated and compared to ex- 
periment. For this we must finally specify the hybridization 
strength T; in the following we take F = 0.5meV (noting that 
comparison to experiment is not critically dependent on this 
choice, with values in the range ~ 0.3 — 0.6meV being found 
quite acceptable). The resultant differential conductance map 
for T = 30mK is shown in fig. [16] (bottom panel), and is in 
rather good agreement with the experimental results reported 
in fig. 1 of fllll. 

In addition to the clear zero-bias Kondo ridges associated 
with both the T and S phases, the conductance map shows 
other features noted in experiment!!. In particular, looking at 
the far left side of the conductance plot, one sees two dark 
'ridges' positioned symmetrically around V s d = 0. As AV g is 
increased the two ridges move together, until they merge to 
form the zero-bias Kondo ridge associated with the T phase. 
The latter persists for a range of AV g , and then the two ridges 
separate again (the pattern being in other words symmetrical 
about AVg. m = — lOmV, for the reasons explained following 
eqn. d62l)). The obvious question arises as to the origin of 
these 'ridges', which we now consider. 

As seen in fig. [16] the T = zero-bias conductance in- 
creases on passing from the S (FL) to the T (USC) phase; and 
in sec. IV Al we showed that such behavior was indicative of 
a vanishing Kondo antiresonance in the single-particle spec- 
trum, as the transition is approached from the S side. This 
is the origin of the ridges seen in the conductance maps, as 
now shown by considering cuts through the conductance map 
of fig. [16] for a sequence of different fixed gate voltages AV 2 . 
The top panel in fig. [T7] accordingly shows the conductance 
vs bias voltage Vsd, for five different values of AV g : — 10m V 
(at the midpoint of the T phase), and AV g = —6, —5, —3 and 
+2 mV as one moves into the S phase (which at T = occurs 
for AVg > — 6.5mV). For AV g = — 10m V the conductance nat- 
urally peaks at V s( j = 0. But on entering the S phase a clear 
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FIG. 17: Top panel: ^-dependence of conductance, taking cuts 
through the conductance map of fig. [16] (bottom) for a sequence 
of different fixed gate voltages: AV g = — lOmV (midpoint of the T 
phase, sold line), and AV g /mV = —6 (long dash), —5 (short dash), 
—3 (point dash) and +2 (dotted), on moving into and through the S 
phase. A clear antiresonance develops in the S phase, as discussed in 
text. The peaks in the conductance, symmetrical about K sc j = 0, lie at 
the center of the ridges seen in the conductance map. 
Bottom panel: Focus on the AV g = — 3mV case. Solid line: conduc- 
tance as in top panel. Long dash: corresponding T = conductance, 
proportional to [D ee (ft) = +jeV s( j) +D ee (co = — ieV s( j)]. Short dash: 
contribution from D ee (co = +^eV s ^) alone, showing a clear Kondo 
antiresonance in the single-particle spectrum itself. 



antiresonance in the conductance is seen to develop, just set- 
ting in by AV g = — 6mV and deepening progressively as AV g 
is increased in the S phase towards — 3mV (then naturally dis- 
appearing as one gets considerably further into the S phase, 
as illustrated by the AV, = +2mV example). And the peaks 
in these conductance profiles, symmetrically disposed about 
V s d = 0, lie at the center of the ridges in the conductance map. 

To show that the above behavior indeed reflects a Kondo 
antiresonance in the single-particle spectrum itself, the bot- 
tom panel in fig. [l7]focusses on the AV, = — 3mV example. 
The solid line again gives the conductance shown in the top 
panel; while the long dash line shows the corresponding T = 
conductance obtained from eqn. doTl ). and thus being propor- 
tional to [D ee {(0 = -fjeVsd) +D ee {(0 = -jcVgd)]. The lat- 
ter clearly captures well the former (the differences naturally 
being due to thermal smearing). Because the conductance is 
proportional to the symmetrized spectra at co = i^eVsd, it is 
not a priori clear that the single-particle spectrum itself con- 
tains a Kondo antiresonance. That it does, however, is seen 
from short-dash line in fig. [T7l bottom), which shows the con- 
tribution from D ee (a> = +\eV s &) itself, seen to contain a clear 
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FIG. 18: As discussed in text, differential conductance map in the 
(V S< 1, AV g )-plane; shown for T = 30mK with the same bare param- 
eters as figs. J16ll7t . except for a slight change in Jn. In this case 
no S/T transition occurs, the system remaining throughout in the S 
phase. The zero-bias Kondo ridge associated with the T phase is thus 
absent, and instead a zero-bias conductance antiresonance persists. 



Kondo antiresonance centered on the Fermi level. We add too 
that, since the peaks/ridges in the conductance stem from the 
'peaks' inherent to the Kondo antiresonance in D ee , they are 
obviously not interpretable in terms of isolated dot states. 

Finally, the zero-bias Kondo ridge in the conductance map 
- formed as described above on merging of the V st j ^ con- 
ductance ridges, and concomitant vanishing of a conductance 
antiresonance (as in fig. [l7]top) - reflects of course the ex- 
istence of the T (USC) phase, and hence a transition to it 
from the S (FL) phase. However one can readily envisage 
a situation where the underlying bare parameters of the sys- 
tem/device are slightly different, such that on ramping down 
the gate voltage the resultant trajectory y — x + Ae comes close 
to but 'misses' the S/T transition; the system as such always 
remaining in a S phase (see e.g. fig. 0a)). In this case no zero- 
bias Kondo ridge associated with the T phase can arise. In- 
stead, from the discussion above, one might intuitively expect 
continued persistence of the conductance antiresonance on de- 
creasing AVo, with attendant finite-bias conductance ridges 
which never quite merge together. That this indeed occurs 
is illustrated in fig. [18] where the conductance map (here for 
T = 30mK) is shown for the same bare parameters as figs. 
( 1161171 ). except for a slight change in /h to —0.6 (the same be- 
havior arising also on changing e.g. Ae rather than /h). And 
the qualitative behavior seen here is indeed similar to that ob- 
served in a second device, shown in fig. 3 of iTTTh (although in 
this case we have not made a quantitative comparison). 



VII. CONCLUDING REMARKS 

As exemplars of multilevel quantum dot systems, we have 
considered in this paper correlated two-level quantum dots, 
coupled in a 1 -channel fashion to metallic leads. Thermo- 
dynamics, single-particle dynamics and electronic transport 
properties show the physical behavior of the system to be 



rich and varied; and our aim has been to obtain a unified 
understanding of the problem for essentially arbitrary dot 
charge/occupancy. Excepting points of high symmetry where 
first order level-crossing transitions arise, associated quantum 
phase transitions are of Kosterlitz-Thouless type, evident in a 
vanishing Kondo scale as the transition to the underscreened 
spin-1 phase is approached from the Fermi liquid side; and 
manifest in particular by a discontinuous jump in the zero-bias 
conductance as the transition is crossed, which we have shown 
can be understood here from an underlying Friedel-Luttinger 
sum rule. We add in fact that an abrupt conductance change 
appears to be a general signature of a KT transition, such be- 
havior arising generally not only in the present model, but also 
in capacitively coupled 2-channel quantum dots which exhibit 
a KT transition from a charge-Kondo Fermi liquid state (with 
a quenched charge pseudospin) to a non-Fermi liquid, doubly 
degenerate 'charge-ordered' phased; and in the problem of 
spinless, capacitively coupled metallic islands/large dots close 
to a degeneracy point between N and N + 1 electron states, de- 
scribed by two Ising-coupled Kondo impurities^. 

Several issues naturally remain to be addressed. We believe 
for example that the generalization of Luttinger's theorem to 
the singular Fermi liquid USC phase (sec. UVBl is significant, 
and raises important basic questions (such as why, and what 
fundamentally does it reflect?). While we do not doubt its 
validity, we have however demonstrated it only numerically; 
and a proper analytical understanding of the result is obvi- 
ously desirable. In this work we have also considered the 
system in the absence of an applied magnetic field, B. Inter- 
esting physics arises also for fi^O (see e.g. where the 
underlying quantum phase transitions are naturally smeared 
into crossovers. In fact, for the USC phase the limits of zero 
field and B — > 0+ are different for T = 0, reflecting the to- 
tal polarization of a free spin- 1/2 (as for the USC fixed point) 
on application of even an infinitesimal field. We will turn in 
subsequent work to the effects of magnetic fields upon single- 
particle dynamics and transport in the model. 
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APPENDIX A: EFFECTIVE LOW-ENERGY MODELS 

We first sketch the derivation of the effective low-energy 
model considered in sec. IIIIB1 spanned by the (1,1) triplet 
and (1, 1) singlet states of the isolated dot; with energies under 
H D of E T = £\ + e 2 + U' - jJh and E s = E\ + £ 2 + U' + |/ H 
respectively. The local unity operator for the dot is 

1 = 1 T + Is (Al) 
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with It = Y*s~ \$ = ^-}S Z )(S = Ij^I m an obvious notation, 
and likewise 1$ = |0,0)(0,0|. These satisfy the following 
identities in the local Hilbert space, 



In 



Si -s 2 



■1 



U 



■ Sl • §2 



(A2) 



as follows using §i ■ §2 = 5 (§ 2 — |) with S = §1 + §2. 

Omitting for brevity the lead contribution Z/l (eqn. ©), the 

low-energy model is given by H e ff = E$ Is +£Y It + "eff ■ The 
first two terms are simply the bare energies of the dot states; 
using eqns. (1A1IA21 they may be written as £5 Is +£Y1t = 
j {E$ + 3E T ) 1 + (E T — E s )si ■ §2, or equivalently as — 7h§i • §2 
on omitting the first (constant/common) term: 



E 



off 



-/ H §1 '§2 + H 



(2) 



off 



(A3) 



Here is the leading (&(V 2 )) contribution arising from 
tunnel coupling to the leads (eqn. (f2]i with V2 = V\ = V, here 
denoted as ft'), given from a SW transformation^ as 



H'l a 



(A4) 

with oc,/3 G {5, r} (retardation effects are as usual neglected). 

In analyzing eqn. (1A4I) one encounters the following 'natu- 
ral' exchange couplings Jf > 0, 



ja = Ny 2 



A£, c 



AE a 



(A5) 



(with N the number of k-states in the lead, such that NV 2 ~ 
<£?(1)). Here the AE" > denote electron removal excitation 
energies from level ; = 1 or 2, relative to the a = T or S dot 
ground state, and the AE" > correspondingly denote elec- 
tron addition energies to level i relative to the a = T or S 
ground state; e.g. AE( — Ed(0, 1) — Ej or AE| = ££>(1,2) — 
£j in the notation of sec. Ill Bl Denoting 



hu — 4/h 



(A6) 



these excitation energies are easily shown to be given by 



AE?=X a -; 



(A7a) 



AE% = l a -y AE 2 a = A«+y 

1 



where (eq n. JTH ) x = £i + jf/ + 1/', y = £2 + W + U'. Hence 
from eqn. 



(A8) 



Jf = J a {x) J? = J a (y) 
with J a (x) defined by: 

1 1 



J a {x) = 



Xa — ■ 



X a +X 



= J a {-x) (A9) 



Direct analysis of eqn. dA4l ) yields, after a standard if labo- 
rious calculation, the effective low-energy model eqn. ( IToT ). 

ft eB = Ji(x,y) §1 -so + J2(x,y) s 2 -s - I fay) §1 -§2 (A 10) 



where potential scattering contributions are omitted for clar- 
ity, and §0 denotes the spin density of the conduction channel 
at the dot. The direct exchange coupling between spins 1 and 
2 is found to be 

I fay) = J H + \ [j T ( X )+J T (y)-J s (x)~J s (y)] (All) 

while the 0{V 2 ) antiferromagnetic exchange couplings be- 
tween spins 1 or 2 and the lead are given by 

hfay) = ±[3J T (x)+J T (y) + (j s (x)~J s (y))] (A12a) 



J 2 fay) = Jifax). (Al2b) 

Using eqn. dA9K these exchange couplings satisfy 

Jifay) = Ji(-x,-y) :i=l,2 (A13a) 



!( x iy) = !(-x,-y) = i(y,x) 



(A13b) 



by virtue of which i7 e ff = Hefifay) m e q n - dAlOb satisfies 
H e s(x,y) = H e g(—x,—y) (reflecting its symmetry under a 
p-h transformation); while from eqn. ( IA1 2bb . ff e ff(^,y) — > 
H e s(y,x) under the 1-2 transformation (eqn. (fT3ll), as expected 
on general grounds from sec. Ill Al These symmetries, and the 
consequent invariance of the phase boundaries in the (x,y)- 
plane to both inversion and reflection about the line y = x, 
are also naturally satisfied when potential scattering terms, 
omitted explicitly from eqn. (|A10t , are included. By contrast, 
the apparent reflection symmetries H e ff(x,y) = H e s(—x,y) = 
H e ff(x, —y) which hold for eqn. ( lAlOb itself (via eqn. JA9b >. 
are not preserved when potential scattering is included; which 
is why the phase boundaries in fig. [7] are not invariant to re- 
flection about the lines x = and v = 0. 

Along the lines y = ±x in the (x,y)-plane, it follows directly 
from eqns. ( |A12b|A13ab that J2 (x, ±x) = J\ (x, ±x), and hence 
from eqns. (|A12a|A9b 



J2{x,±x) = J\(x,±x) = 2J T (x), 



(A14) 



with I(x, ±x) = 7h + [J fa) — J S (x)} following similarly from 
eqn. dAl lb . In consequence, as discussed in sec. HIIBl i/ e ff is 
separable along the lines y = ±x, and first order level-crossing 
transitions thus arise. Note incidentally that for given 7h < 
(A7b) 0' b° tri J T ( X ) an d I(x,±x) increase on increasing |x| from 0, 



both of which act (see eqn. ( fTTI i) to favor the triplet phase; 
consistent with the USC phase surviving for \x\ > after it 
has been destroyed for x = 0, as indeed found in fig. [7] 

Finally, a byproduct of the above gives directly the pure 
spin-1 Kondo model appropriate for ferromagnetic /h > 
deep in the USC spin-1 phas o 15 ' 16 , where the (1,1) singlet 
dot states are energetically irrelevant and only the (1,1) triplet 
states need be retained. To this end simply project into the 
triplet sector, H K = lT^efflT- Writing 7i§i • §0 +/2S2 ■ §0 = 
\{J\ +h){$\ +§2) -s +\{J\ -h){%\ -82) -8o. recognising 
that 1t(§i — §2) It = and neglecting constants, eqn. dAl 01 ) 
gives a spin-1 Kondo model 



H K = l T HefflT = /*rS-§o 



(A15) 
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where S = 1t(§i +§2)1t is a P ure spin-1 operator, and Jk = 
j(Ji +J2) is given from eqn. ( IA125 by 

J K = J T {x)+J T (y) (A16) 

with J T (x) given explicitly by eqn. dA9t , Along the lines 
y = ±x in particular, eqns. (IA16IA9IA6I > give Jk = 2J T (x); 
i.e. pJ K = fr[t/ + \Jyi]/{[U + ±/ H ] 2 -x 2 ) with p here the 
lead density of states per conduction orbital (such that F = 
nV 2 Np). From perturbative scaling 1 the spin-1 Kondo scale 
Tg follows as 7jf ~ Dexp(— 1/pJk) (with the exponen- 
tial dependence as usual of the essence, and the prefactor im- 
material), and eqn. d42b for T^ =i thus results. From NRG 
calculations we have confirmed explicitly that the dependence 
of 7Jf =1 on U + \ /h is indeed as predicted by eqn. (142) . 

APPENDIX B: SELF-ENERGIES 

The key NRG method for calculating the self-energy 4 -^ is 
readily extended to multilevel dots/impurities. With H\ the 
interacting part of the dot Hamiltonian, equation of motion 
techniques^ 2 - are used to obtain the following basic equation 
for the retarded propagators {Gjj((o)}: 

£((©+ ■-e i )8 u -r a (<o))G,j(<o) = Sij + dld^HjU^}) 

(Bl) 

The sum is over the dot levels (/ = 1,2 here), and [, ] denotes 
a commutator. By definition, Hj = in the non-interacting 
limit; whence eqn. (IB 11 1 is of form 

[G°(fl>)] _1 G(a>) = 1 + F(fi>) (B2) 

with G°((o) the non-interacting propagator matrix and the el- 
ements of F(co) given by: 

Ftj(a>) = ((feAMja)) < B3 > 

Using the Dyson equation in the form [G ] -1 = [G] _1 + S, 
eqn. ( IB2b gives directly eqn. d39l l, 

S(ffl) = F{co)[G(co)}- 1 (B4) 



from which the self-energies are calculated directly (sec. [V]). 
Combining eqn. dB4) with the Dyson equation also gives G = 
G°(l + F),sothat eqn. (IB4I) may be written alternatively as 

S(o) = F(co) [l + Fico)}- 1 [G°(co)Y l (B5) 

(which we exploit to calculate the renormalized level e* when 
considering dynamics on the y = x line, eqn. ( 1541 sec. IV B| ). 

Hi is given explicitly for the present problem by the sep- 
arable sum (eqn. (Q~|)) Hj = U^A^n^ + U'h x h 2 — Jh §1 • §2. 
The elements eqn. (IB3b of F are thus linearly separable as 
Fjj = Fjj +F% +Fji (in obvious notation), and are calculated 
individually. Since each such term is a retarded correlation 
function, they are Lehmann resolvable^! and in consequence 
satisfy sum rules. Writing Fij(co) = F/j((o) —iF^co) (with 
the real/imaginary parts related by Hilbert transformation), the 
general sum rule is 

J^Ffj(co) = {{[d ic ,HjU} a }) (B6) 

where {,} denotes an anticommutator. Specifically, for the 
present problem in the absence of an applied magnetic field, 
it is readily shown that the diagonal-element sum rule is 

/Lt f * (<d) = * utfti) + ^ (B7) 

(where i means the opposite level to /); while for the off- 
diagonal elements: 

jT ^ Ffj(co) & (-W + |/ H ) {JjM 

= (-U'+pn) J° d(0Dij((0) (B8b) 

These sum rules provide a check on the accuracy of the NRG 
calculations, and in practice are well satisfied. 
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